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I. INTRODUCTION AND OUTLINE 


Most aspects of high energy physics beyond the standard model can only be tested 
by going to energies far greater than those which present accelerators can provide. 
Fortunately, the marriage between particle physics and cosmology has provided a way 
to “experimentally” test the new theories of fundamental forces. 

The Iffiy realization, discovered both in the context of the inflationary Universe 
scenario B and of topological defects models 0 is that physics of the very early Uni¬ 
verse may explain the origin of structure in the Universe. It now appears that a 
rich set of data concerning the nonrandom distribution of matter on a wide range 
of cosmological scales, and on the anisotropies in the cosmic microwave background 
(CMB), may potentially be explained by high energy physics. In addition, studying 
the consequences of particle physics models in the context of cosmology may lead 
to severe constraints on new microscopic theories. Finally, particle physics and held 
theory may provide explanations of some deep cosmological puzzles, e.g. why the 
Universe at the present time appears so homogeneous, so close to being spatially hat, 
and why it contains the observed small net baryon to entropy ratio. 

In these lectures, I focus on three important aspects of modern cosmology. The 
hrst concerns some fundamental problems of inhationary cosmology. In particular, 
some recent progress in the understanding of “reheating” in inhation will be reviewed. 

The second topic is the classical and quantum theory of cosmological perturba¬ 
tions, the main tool of modern cosmology. A general relativistic and quantum me¬ 
chanical analysis of the generation and evolution of linearized huctuations is essential 
in order to be able to accurately calculate the amplitude of density perturbations and 
CMB anisotropies. 

As a third topic, I discuss the role of topological defects in baryogenesis and as 
possible sees for cosmological structure formation. 

The specihc outline is as follows: 

1. Introduction and Outline 

2. Lecture 1: Inflationary Universe: Progress and Problems 

2.A Problems of Standard Cosmology 
2.B Inflationary Universe Scenario 
2.C Problems of Inflation 
2.D Inflation and Nonsingular Cosmology 
2.E Reheating in Inflationary Cosmology 

2. F Summary 

3. Lecture 2: Classical and Quantum Theory of Cosmological Perturba¬ 
tions 

3. A Basic Issues 

3.B Newtonian Theory 
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3.C Relativistic Theory: Classical Analysis 
3.D Relativistic Theory: Quantum Analysis 

3. E Summary 

4. Lecture 3: Topological Defects, Structure Formation and Baryogen- 
esis 

4. A Quantifying Data on Large-Scale Structure 
4.B Topological Defects 

4.C Formation of Defects in Cosmological Phase Transitions 

4.D Evolution of Strings and Scaling 

4.E Cosmic Strings and Structure Eormation 

4.E Specihc Predictions 

4.G Principles of Baryogenesis 

4.H GUT Baryogenesis and Topological Defects 

4.1 Electroweak Baryogenesis and Topological Defects 

4.J Summary 

Unless otherwise specihed, units in which h = c = ks = ^ will be used. Distances 
are expressed in Mpc (Ipc ~ 3.06 light years). Eollowing the usual convention, h 
indicates the expansion rate of the Universe in units of 100 km s“^ Mpc“^, D = p/Pc 
is the ratio of the energy density p to the critical density pc (the density which yields 
a spatially flat Universe), G is Newton’s constant and rupi is the Planck mass. 


II. INFLATIONARY UNIVERSE: PROGRESS AND PROBLEMS 

The hypothesis that the Universe underwent a period of exponential expansion at 
very early times has become the most popular theory of the early Universe. Not only 
does it solve some of the problems of standard big bang cosmology, but it also provides 
a causal theory for the origin of inhomogeneities in the Universe which is predictive 
and in reasonable agreement with current observational results. Nevertheless, there 
are several problems of principle which merit further study. 


A. Problems of Standard Cosmology 


The standard big bang cosmology rests on three theoretical pillars: the cosmolog¬ 
ical principle, Einstein’s general theory of relativity and a perfect fluid description of 
matter. 

The cosmological principle states that on large distance scales the Universe is 
homogeneous. This implies that the metric of space-time can be written in Eriedmann- 
Robertson-Walker (ERW) form: 


ds^ 


a{tY 


dr'^ 

1 — kr‘^ 


+ r‘^{d'd‘^ + sin^ iddp'^) , 


( 1 ) 
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where the constant k determines the topology of the spatial sections. In the following, 
we shall usnally set k = 0, i.e. consider a spatially closed Universe. In this case, we 
can withont loss of generality take the scale factor a(t) to be eqnal to 1 at the present 
time to, i.e. a(to) = 1. The coordinates r, {} and ip are comoving spherical coordinates. 
World lines with constant comoving coordinates are geodesics corresponding to par¬ 
ticles at rest. If the Universe is expanding, i.e. a{t) is increasing, then the physical 
distance Axp{t) between two points at rest with hxed comoving distance Axc grows: 

Axp = a{t)Axc . ( 2 ) 

The dynamics of an expanding Universe is determined by the Einstein eqnations, 
which relate the expansion rate to the matter content, specihcally to the energy 
density p and pressnre p. For a homogeneons and isotropic Universe, they rednce to 
the Friedmann-Robertston-Walker (FRW) eqnations 

= ( 3 ) 

\a/ 6 


d AttG , ^ 

- = -^(P + 3p). (4) 

a 6 

These eqnations can be combined to yield the continnity eqnation (with Hnbble con¬ 
stant H = a/a) 


p = -3H{p + p). 


( 5 ) 


The third key assnmption of standard cosmology is that matter is described by 
an ideal gas with an eqnation of state 


p = wp. (6) 

For cold matter, pressnre is negligible and hence ta = 0. From (^) it follows that 

Pm{t) ~ a"^(t), (7) 

where pm is the energy density in cold matter. For radiation we have w = 1/3 and 
hence it follows from that 

Pr{t) ^ , ( 8 ) 

Pr{t) being the energy density in radiation. 

The three classic observational pillars of standard cosmology are Hnbble’s law, the 
existence and black body natnre of the nearly isotropic CMB, and the abnndances 
of light elements (nncleosynthesis). These snccesses are discnssed in detail in many 
textbooks on cosmology, and will therefore not be reviewed here. 
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It is, however, important to recall two important aspects concerning the thermal 
history of the early Universe. Since the energy density in radiation redshifts faster 
than the matter energy density, it follows by working backwards in time from the 
present data that although the energy density of the Universe is now mostly in cold 
matter, it was initially dominated by radiation. The transition occurred at a time 
denoted by teq, the “time of equal matter and radiation”. As will be discussed in 
Section 3, teq is the time when perturbations can start to grow by gravitational 
clustering. The second important time is Ueo the “time of recombination” when 
photons fell out of equilibrium (since ions and electrons had by then combined to 
form electrically neutral atoms). The photons of the CMB have travelled without 
scattering from tree- Their spatial distribution is predicted to be a black body since the 
cosmological redshift preserves the black body nature of the initial spectrum (simply 
redshifting the temperature) which was in turn determined by thermal equilibrium. 
CMB anisotropies probe the density fluctuations at t^ec- Note that for the usual 
values of the cosmological parameters, teq < tree- 

Standard Big Bang cosmology is faced with several important problems. Only one 
of these, the age problem, is a potential conflict with observations. The others which 
I will focus on here - the homogeneity, flatness and formation of structure problems 
(see e.g. [Q) - are questions which have no answer within the theory and are therefore 
the main motivation for the new cosmological models which will be discussed in later 
sections of these lecture notes. 

The horizon problem is illustrated in Fig. 1. As is sketched, the comoving region 
^p{trec) over which the CMB is observed to be homogeneous to better than one part in 
10^ is much larger than the comoving forward light cone if{trec) at tree, which is the 
maximal distance over which microphysical forces could have caused the homogeneity: 




( 10 ) 


0 


From the above equations it is obvious that ip{trec) 3> ifiUec)- Hence, standard 
cosmology cannot explain the observed isotropy of the CMB. 


5 


1 

i 



FIG. 1. A space-time diagram (physical distance Xp versus time t) illustrating the 
homogeneity problem: the past light cone ip{t) at the time tree of last scattering is much 
larger than the forward light cone if{t) at tree- 


In standard cosmology and in an expanding Universe, = 1 is an nnstable fixed 
point. This can be seen as follows. For a spatially flat Universe (U = 1) 




SttG 

7^ Pc ) 


whereas for a nonflat Universe 


H^ + eT^ 


SttG 



(aTf ■ 


( 11 ) 


( 12 ) 


(13) 


The qnantity £ is proportional to s where s is the entropy c 
standard cosmology, e is constant. Combining ([TT|) and (|T^ gives 


Hence, in 


p- Pc _ 3 sT"^ 

Pe SnG pe 


(14) 


Thns, as the temperatnre decreases, H — 1 increases. In fact, in order to explain the 
present small valne of H — 1 ~ 0(1), the initial energy density had to be extremely 
close to critical density. For example, at T = 10^® GeV, (JT^) implies 


pc 

What is the origin of these hne tnned initial conditions? This is the flatness problem 
of standard cosmology. 
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The third of the classic problems of standard cosmological model is the “forma¬ 
tion of strnctnre problem.” Observations indicate that galaxies and even clusters of 
galaxies have nonrandom correlations on scales larger than 50 Mpc (see e.g. p,^). 
This scale is comparable to the comoving horizon at teg. Thus, if the initial density 
perturbations were produced much before teq, the correlations cannot be explained by 
a causal mechanism. Gravity alone is, in general, too weak to build up correlations 
on the scale of clusters after teq (see, however, the explosion scenario of |^). Hence, 
the two questions of what generates the primordial density perturbations and what 
causes the observed correlations, do not have an answer in the context of standard 
cosmology. This problem is illustrated by Fig. 2. 


t 



FIG. 2. A sketch (conformal separation vs. time) of the formation of structure problem: 
the comoving separation dc between two clusters is larger than the forward light cone at 
time teq- 


There are other serious problems of standard cosmology, e.g. the age and the 
cosmological constant problems. However, to date modern cosmology does not shed 
any light on these problems, and I will therefore not address them here. 


B. Inflationary Universe Scenario 

The idea of inflation 0 is very simple (for some early reviews of inflation see e.g. 
P-0). We assume there is a time interval beginning at ti and ending at tji (the 
“reheating time”) during which the Universe is exponentially expanding, i.e., 

a{t) ~ te [U, tn] (16) 

with constant Hubble expansion parameter H. Such a period is called “de Sitter” 
or “inflationary.” The success of Big Bang nucleosynthesis sets an upper limit to the 
time of reheating: 

tn <C tj^s , (17) 

tNs being the time of nucleosynthesis. 
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FIG. 3. The phases of an inflationary Universe. The times U and tn denote the beginning 
and end of inflation, respectively. In some models of inflation, there is no initial radiation 
dominated FRW period. Rather, the classical space-time emerges directly in an inflationary 
state from some initial quantum gravity state. 

The phases of an inflationary Universe are sketched in Fig. 3. Before the onset 
of inflation there are no constraints on the state of the Universe. In some models a 
classical space-time emerges immediately in an inflationary state, in others there is 
an initial radiation dominated FRW period. Our sketch applies to the second case. 
After tR, the Universe is very hot and dense, and the subsequent evolution is as 
in standard cosmology. During the inflationary phase, the number density of any 
particles initially in thermal equilibrium at t decays exponentially. Hence, the 
matter temperature Tmit) also decays exponentially. At t = tR, all of the energy 
which is responsible for inflation (see later) is released as thermal energy. This is a 
nonadiabatic process during which the entropy increases by a large factor. 

Fig. 4 is a sketch of how a period of inflation can solve the homogeneity problem. 
At = tR—ti is the period of inflation. During inflation, the forward light cone increases 
exponentially compared to a model without inflation, whereas the past light cone is 
not affected for t > tR. Hence, provided At is sufficiently large, if(tR) will be greater 
than ipitR). 
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FIG. 4. Sketch (physical coordinates vs. time) of the solution of the homogeneity 
problem. During inflation, the forward light cone lf{t) is expanded exponentially when 
measured in physical coordinates. Hence, it does not require many e-foldings of inflation in 
order that Ij (t) becomes larger than the past light cone at the time of last scattering. The 
dashed line is the forward light cone without inflation. 


Inflation also can solve the flatness problem MEJ The key point is that the entropy 
density s is no longer constant. As will be explained later, the temperatures at 
and tR are essentially equal. Hence, the entropy increases during inflation by a factor 
exp(3ifAt). Thus, e decreases by a factor of exp{—2HAt). Hence, (p — Pc)/p can 
be of order 1 both at and at the present time. In fact, if inflation occurs at all, 
then rather generically, the theory predicts that at the present time H = 1 to a l^h 
accuracy (now H < 1 requires special initial conditions or rather special models H). 

Most importantly, inflation provides a mechanism which in a causal way generates 
the primordial perturbations required for galaxies, clusters and even larger objects. 
In inflationary Universe models, the Hubble radius (“apparent” horizon), 3t, and the 
“actual” horizon (the forward light cone) do not coincide at late times. Provided that 
the duration of inflation is sufficiently long, then (as sketched in Fig. 5) all scales 
within our apparent horizon were inside the actual horizon since tj. Thus, it is in 
principle possible to have a casual generation mechanism for perturbations HhUj. 

The generation of perturbations is supposed to be due to a causal microphysical 
process. Such processes can only act coherently on length scales smaller than the 
Hubble radius init) where 




( 18 ) 


A heuristic way to understand the meaning of init) is to realize that it is the distance 
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which light (and hence the maximal distance any causal effects) can propagate in one 
expansion time. 


f 



FIG. 5. A sketch (physical coordinates vs. time) of the solution of the formation of 
structure problem. Provided that the period of inflation is sufficiently long, the separation 
dc between two galaxy clusters is at all times smaller than the forward light cone. The 
dashed line indicates the Hubble radius. Note that dc starts out smaller than the Hubble 
radius, crosses it during the de Sitter period, and then reenters it at late times. 


As will be discussed in Chapter 4, the density perturbations produced during 
inflation are due to quantum fluctuations in the matter and gravitational fields ollJ . 
The amplitude of these inhomogeneities corresponds to a temperature Th 

Th ^ H , (19) 

the Hawking temperature of the de Sitter phase. This implies that at all times t during 
inflation, perturbations with a fixed physical wavelength ~ H~^ will be produced. 
Subsequently, the length of the waves is stretched with the expansion of space, and 
soon becomes larger than the Hubble radius. The phases of the inhomogeneities 
are random. Thus, the inflationary Universe scenario predicts perturbations on all 
scales ranging from the comoving Hubble radius at the beginning of inflation to 
the corresponding quantity at the time of reheating. In particular, provided that 
inflation lasts sufficiently long, perturbations on scales of galaxies and beyond will be 
generated. Note, however, that it is very dangerous to interpret de Sitter Hawking 
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radiation as thermal radiation. In fact, the equation of state of this “radiation” is 
not thermal 0. 

Obviously, the key question is how to obtain inflation. From the FRW equations, 
it follows that in order to get exponential increase of the scale factor, the equation of 
state of matter must be 


P= -P (20) 

This is where the connection with particle physics comes in. The energy density and 
pressure of a scalar quantum held tp are given by 

pM = + + VM ( 21 ) 

pM = \p - - v(<fi). ( 22 ) 

Thus, provided that at some initial time ti 

ti) = V(p{Xiti) = 0 (23) 

and 

V{p{xi, U)) > 0 , (24) 


the equation of state of matter will be (|^) . 

The next question is how to realize the required initial conditions 
maintain the key constraints 


and to 


« V(,p) , (VvY « V(,p) 


(25) 


for sufficiently long. Various ways of realizing these conditions were put forward, and 
they gave rise to different models of inhation. I will focus on “old inhation,” “new 
inhation”” and “chaotic inhation.” There are many other attempts at producing an 
inhationary scenario, but there is as of now no convincing realization. 

Old Inflation 

The old inflationary Universe model is based on a scalar held theory which 
undergoes a flrst order phase transition. As a toy model, consider a scalar held 
theory with the potential V{(p) of Figure 6. This potential has a metastable “false” 
vacuum at (p = 0, whereas the lowest energy state (the “true” vacuum) is (p = a. 
Finite temperature ehects 0 lead to extra terms in the flnite temperature ehective 
potential which are proportional to (the resulting flnite temperature effective 

potential is also depicted in Figure 6). Thus, at high temperatures, the energetically 
preferred state is the false vacuum state. Note that this is only true if p is in thermal 
equilibrium with the other fields in the system. 
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FIG. 6. The finite temperature effective potential in a theory with a first order phase 
transition. 


For fairly general initial conditions, (p{x) is trapped in the metastable state ip = 0 
as the Universe cools below the critical temperature Tc. As the Universe expands 
further, all contributions to the energy-momentum tensor except for the contri¬ 
bution 


(26) 

redshift. Hence, provided that the potential V(p) is shifted upwards such that V (a) = 
0, then the equation of state in the false vacuum approaches p = —p, and inflation 
sets in. Aft^ a period F”^, where F is the tunnelling rate, bubbles of <^9 = a begin 
to nucleate Ej in a sea of false vacuum p = 0. Inflation lasts until the false vacuum 
decays. During inflation, the Hubble constant is given by 

w" = ^ 1/(0). (27) 

Note that the condition V (a) = 0, which looks rather unnatural, is required to avoid 
a large cosmological constant today (none of the present inflationary Universe models 
manages to circumvent or solve the cosmological constant problem). 

It was immediately realized that old inflation has a serious “graceful exit” 
problem 00. The bubbles nucleate after inflation with radius r -C and would 
today be much smaller than our apparent horizon. Thus, unless bubbles percolate, 
the model predicts extremely large inhomogeneities inside the Hubble radius, in con¬ 
tradiction with the observed isotropy of the microwave background radiation. 

For bubbles to percolate, a sufficiently large number must be produced so that 
they collide and homogenize over a scale larger than the present Hubble radius. How¬ 
ever, with exponential expansion, the volume between bubbles expands exponentially 
whereas the volume inside bubbles expands only with a low power. This prevents 
percolation. 

New Inflation 
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Because of the graceful exit problem, old inflation never was considered to be a 
viable cosmological model. However, soon after the seminal paper by Guth, Linde 0 
and independently Albrecht and Steinhardt il put forwards a modihed scenario, the 
New Inflationary Universe. 

The starting point is a scalar held theory with a double well potential which 
undergoes a second order phase transition (Fig. 7). V{(p) is symmetric and = 0 is 
a local maximum of the zero temperature potential. Once again, it was argued that 
hnite temperature effects conhne (p{x) to values near (p = 0 at temperatures T > T^. 
For T < Tc, thermal huctuations trigger the instability of (p{x) = 0 and (p{x) evolves 
towards either of the global minima at = ±<7 by the classical equation of motion 

(p + 3Hlp — a~‘^ ‘■P = ■ (28) 


Within a huctuation region, (p{x) will be homogeneous. In such 
neglect the spatial gradient terms in Eq. (^81). Then, from (plj) and 
off the induced equation of state. The condition for inhation is 


a region, we can 
(^^ we can read 


•p" « V(p), 


(29) 


i.e. slow rolling. Often, the “slow rolling” approximation is made to hnd solutions 
of (pSl). This consists of dropping the p term. 



FIG. 7. The finite temperature effective potential in a theory with a second order phase 
transition. 


There is no graceful exit problem in the new inflationary Universe. Since the 
fluctuation domains are established before the onset of inflation, any boundary walls 
will be inflated outside the present Hubble radius. 

Let us, for the moment, return to the general features of the new inflationary 
Universe scenario. At the time tc of the phase transition, p{t) will start to move 
from near (p = 0 towards either ±ct as described by the classical equation of motion, 
i.e. (^8[) . At or soon after tc, the energy-momentum tensor of the Universe will start 
to be dominated by V(p), and inflation will commence, ti shall denote the time of the 
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onset of inflation. Eventnally, 0(t) will reach large valnes for which nonlinear effects 
become important. The time at which this occnrs is ts- For t > ts, (p{t) rapidly 
accelerates, reaches ±cr, overshoots and starts oscillating abont the global minimum 
of V{ip). The amplitnde of this oscillation is damped by the expansion of the Universe 
and (predominantly) by the conpling of p to other helds. At time tn, the energy in p 
drops below the energy of the thermal bath of particles prodnced dnring the period 
of oscillation. 

The evolntion of p{t) is sketched in Fig. 8. The time period between ts and tji 
is called the reheating period and is nsnally short compared to the Hnbble expansion 
time. For t > tji, the Universe is again radiation dominated. 


f T 




FIG. 8. Evolution of p{t) and T{t) in the new inflationary Universe. 


In order to obtain inflation, the potential V{p) mnst be very flat near the false 
vacnnm at <^9 = 0. This can only be the case if all of the conpling constants appearing 
in the potential are small. However, this implies that the p cannot be in thermal 
eqnilibrinm at early times, which wonld be reqnired to localize p in the false vacnnm. 
In the absence of thermal eqnilibrinm, the initial conditions for p are only constrained 
by reqniring that the total energy density in p not exceed the total energy density 
of the Universe. Most of the phase space of these initial conditions lies at valnes of 
\p\ » a. This leads to the “chaotic” inflation scenario il. 

Chaotic Inflation 

Consider a region in space where at the initial time p{x) is very large, homo- 
geneons and static. In this case, the energy-momentnm tensor will be immediately 
dominated by the large potential energy term and indnce an eqnation of state p ~ —p 
which leads to inflation. Dne to the large Hnbble damping term in the scalar held 
eqnation of motion, p{x) will only roll very slowly towards p = The kinetic energy 
contribntion to will remain small, the spatial gradient contribntion will be expo¬ 
nentially snppressed dne to the expansion of the Universe, and thns inflation persists. 
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Note that in contrast to old and new inflation, no initial thermal bath is required. 
Note also that the precise form of V(cp) is irrelevant to the mechanism. In particular, 
V((p) need not be a double well potential. This is a significant advantage, since for 
scalar fields other than Higgs fields used for spontaneous symmetry breaking, there 
is no particle physics motivation for assuming a double well potential, and since the 
inflaton (the held which gives rise to inflation) cannot be a conventional Higgs field 
due to the severe fine tuning constraints. 

The held and temperature evolution in a chaotic inflation model is similar to what 
is depicted in Figure 8, except that (p is rolling towards the true vacuum at (p = a 
from the direction of large held values. 

Chaotic inflation is a much more radical departure from standard cosmology than 
old and new inflation. In the latter, the inflationary phase can be viewed as a short 
phase of exponential expansion bounded at both ends by phases of radiation dom¬ 
ination. In chaotic inflation, a piece of the Universe emerges with an inflationary 
equation of state immediately after the quantum gravity (or string) epoch. 

The chaotic inflationary Universe scenario has been developed in great detail (see 
e.g. |2M for a recent review). One important addition is the inclusion of stochastic 
noise in the equation of motion for (p in order to take into account the effects of 
quantum fluctuations. It can in fact be shown that for sufficiently large values of 
1931 , the stochastic force terms are more important than the classical relaxation force 
V'((p). There is equal probability for the quantum fluctuations to lead to an increase 
or decrease of \p>\. Hence, in a substantial fraction of comoving volume, the held 
(p will climb up the potential. This leads to the conclusion that chaotic inflation is 
eternal. At all times, a large fraction of the physical space will be inflating. Another 
consequence of including stochastic terms is that on large scales (much larger than 
the present Hubble radius), the Universe will look extremely inhomogeneous. 


C. Problems of Inflationary Cosmology 

In spite of its great success at resolving some of the problems of standard cos¬ 
mology and of providing a causal, predictive theory of structure formation, there are 
several important unresolved conceptual problems in inflationary cosmology. I will 
focus on three of these problems, the cosmological constant mystery, the fluctuation 
problem, and the dynamics of reheating. 

Cosmological Constant Problem 

Since the cosmological constant acts as an effective energy density, its value is 
bounded from above by the present energy density of the Universe. In Planck units, 
the constraint on the effective cosmological constant Ae// is (see e.g. [p6| ) 

^ < 10-^22 
ml, 


(30) 
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This constraint applies both to the bare cosmological constant and to any matter 
contribntion which acts as an effective cosmological constant. 

The trne vacnnm valne of the potential V{^ip) acts as an effective cosmological 
constant. Its valne is not constrained by any particle physics reqnirements (in the 
absence of special symmetries). The cosmological constant problem is thns even more 
accnte in inflationary cosmology than it nsnally is. The same nnknown mechanism 
which mnst act to shift the potential (see Fignre 6) snch that inflation occnrs in the 
false vacnnm mnst also adjnst the potential to vanish in the trne vacnnm. 

Snpersymmetric theories may provide a resolntion of this problem, since nnbro- 
ken snpersymmetry forces V{ip) = 0 in the snpersymmetric vacnnm. However, sn- 
persymmetry breaking will indnce a nonvanishing V{ip) in the trne vacnnm after 
snpersymmetry breaking. 

We may therefore be forced to look for realizations of inflation which do not make 
nse of scalar helds. There are several possibilities. It is possible to obtain inflation in 
higher derivative gravity theories. In fact, the hrst model with exponential expansion 
of the Universe was obtained 0 in an gravity theory. The extra degrees of freedom 
associated with the higher derivative terms act as scalar helds with a potential which 
antomatically vanishes in the trne vacnnm. For some recent work on higher derivative 
gravity inhat ion see also . 

Another way to obtain inhation is by making nse of condensates (see 11 and 


30| for diherent approaches to this problem). An additional motivation for following 


this ronte to inhation is that the symmetry breaking mechanisms observed in natnre 
(in condensed matter systems) are indnced by the formation of condensates such as 
Cooper pairs. Again, in a model of condensates there is no freedom to add a constant 
to the ehective potential. 

The main problem of studying the possibility of obtaining inhation using conden¬ 
sates is that the quantum ehects which determine the theory are highly nonpertur- 
bative. In particular, the ehective potential written in terms of a condensate {(p) 
does not correspond to a renormalizable theory and will in general 0 contain terms 
of arbitrary power in {ip). However (see |^), one may make progress by assuming 
certain general properties of the ehective potential. 

Let us 0 consider a theory in which at some time ti a condensate {ip) forms, i.e. 
{ip) = 0 for t < ti and {ip) ^ 0 for t > ti. The expectation value of the Hamiltonian 
H written in terms of the condensate {ip) contains terms of arbitrary powers of {ip): 


(H) = i:(-l)”n!a„(v>)’ 


(31) 


We summarize our ignorance of the nonperturbative physics in the assumption that 
the resulting series is asymptotic, and in particular Borel summable, with coefficients 
a„ oc 1. In this case, we can resum the series to obtain H 


(H) = / 
Jo 


f{t)dt 


,-l/i 


0 titnipi {ip)) 


(32) 
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( 33 ) 


where the function f{t) is related to the coefficients a„ via 

1 / 
a„ = — / 
nl Jo 

The expectation value of the Hamiltonian (H) can be interpreted as the effective 
potential 14// of this theory. The question is under which conditions this potential 
gives rise to inflation. If we regard (ip) as a classical held (i.e. neglect the ultraviolet 
and infrared divergences of the theory), then the dynamics of the model can be read 
off directly from (|3^) , with initial conditions for {(p) at the time close to (ip) = 0. 
It is easy to check that rather generically, the conditions required to have slow rolling 
of (p, namely 

Viripi « V4:87 iV (34) 

V"mli « 247rl/, (35) 

are satished. However, since the potential decays only slowly at large values of {(p) 
and since there is no true vacuum state at hnite values of (ip), the slow rolling con¬ 
ditions are satished for all times. In this case, inhation would never end - an obvious 
cosmological disaster. 

However, (ip) is not a classical scalar held but the expectation value of a condensate 
operator. Thus, we have to worry about diverging contributions to this expectation 
value. In particular, in a theory with symmetry breaking there will often be massless 
excitations which will give rise to infrared divergences. It is necessary to introduce 
an infrared cutoh energy e whose value is determined in the context of cosmology 
by the Hubble expansion rate. Note in particular that this cutoh is time-dependent. 
Ehectively, we thus have a theory of two scalar helds (ip) and e. In this case, the hrst 
of the slow rolling conditions becomes (if £ is expressed in Planck units) 

i’^rripi + ip'^ << 2V . (36) 

The infrared cutoh changes the form of the ehective potential. We assume that 
this change can be modelled by replacing (ip) by {(p)/e. If we (following |^) take 
the infrared cutoh to be 

e{t) = ^^[1 - a{Hm , (37) 

rripi 

where 0 < a << 1 and p is an integer and the time at the beginning of the rolling has 
been set to t = 0, then it can be shown @ that an period of inhation with a graceful 
exit is realized. After the condensate {(p) starts rolling at (ip) ~ 0, inhation will 
commence. As inhation proceeds, e{t) will slowly grow and will eventually dominate 
the energy functional, signaling an end of the inhationary period. From (|37|) it follows 
that inhation lasts until = 1. 
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This analysis demonstrates that it is in principle possible to obtain inflation from 
condensates. However, the model mnst be stndied in mnch more detail before we can 
determine whether it gives a realization of inflation which is free of problems. 

Fluctuation Problem 


A generic problem for all realizations of inflation stndied np to now concerns the 
amplitnde of the density pertnrbations which are induced by quantum fluctuations 
during the period of exponential expansion. From the amplitude of CMB anisotropies 
measured by COBE, and from the present amplitude of density inhomogeneities on 
scales of clusters of galaxies, it follows that the amplitude of the mass fluctuations 
6M/M on a length scale given by the comoving wavenumber k at the time tnik) 
when that scale crosses the Hubble radius in the FRW period is 

^{k,tH{k)) oc 10"^ (38) 

The generation and evolution of fluctuations will be discussed in detail in Section 
3. The perturbations arise during inflation as quantum excitations. Their amplitude 
at the time ti{k) when the scale k leaves the Hubble radius during inflation is given 
by 


6M 

W 


{k,ti{k)) 




V'S^ 

~ In(fe)) 


(39) 


where 6ip is given by the amplitude of the quantum fluctuation of 6ip{k) (note that 
this is a momentum space quantity). While the scale k is outside of the Hubble 
radius, the fluctuation amplitude grows by general relativistic gravitational effects. 
The amplitudes at ti{k) and tuik) are related by 


5M 

W 


{K,tH{k)) 




1 

l+p/p 


6M 


{k,ti{k)) 


(40) 


(see e.g. ^^). Combining (|35|) and (13) and working out the result for the potential 

V{ip) = (41) 

we obtain the result Hizl 


6M 




(42) 


Thus, in order to agree with the observed value (pSD , the coupling constant A must 
be extremely small: 


A < 10"^^ 


(43) 


It has been shown in that the above conclusion is generic, at least for models in 
which inflation is driven by a scalar field. In order that inflation does not produce a too 
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large amplitude of the spectrum of perturbations, a dimensionless number appearing 
in the potential must be set to a very small value. A possible resolution of this 
problem will be mentioned in the following subsection. 

Reheating Problem 

A question which has recently received a lot of attention and will be discussed in 
greater detail in one of the following subsections is the issue of reheating in inflationary 
cosmology. The question concerns the energy transfer between the inflaton and matter 
helds which is supposed to take place at the end of inflation (see Fig. 8). 

According to either new inflation or chaotic inflation, the dynamics of the inflaton 
leads hrst to a transfer of energy from potential energy of the inflaton to kinetic 
energy. After the period of slow rolling, the inflaton tp begins to oscillate about the 
true minimum of V{ip). Quantum mechanically, the state of homogeneous oscillation 
corresponds to a coherent state. Any coupling of p to other helds (and even self 
coupling terms of p) will lead to a decay of this state. This corresponds to the particle 
production. The produced particles will be relativistic, and thus at the conclusion of 
the reheating period a radiation dominated Universe will emerge. 

The key questions are by what mechanism and how fast the decay of the coherent 
state takes place. It is important to determine the temperature of the produced par¬ 
ticles at the end of the reheating period. The answers are relevant to many important 
questions regarding the post-inhationary evolution. For example, it is important to 
know whether the temperature after reheating is high enough to allow GUT baryogen- 
esis and the production of GUT-scale topological defects. In supersymmetric models, 
the answer determines the predicted abundance of gravitinos and other moduli helds. 

Recently, there has been a complete change in our understanding of reheating. 
This topic will be discussed in detail below. 


D. Inflation and Nonsingnlar Cosmology 

The question we wish to address in this subsection is whether it is possible to 
construct a class of ehective actions for gravity which have improved singularity prop¬ 
erties and which predict inflation, with the constraint that they give the correct low 
curvature limit. Since Planck scale physics will generate corrections to the Einstein 
action, it is quite reasonable to consider higher derivative gravity models. 

What follows is a summary of recent work il in which we have constructed an 
ehective action for gravity in which all solutions with sufficient symmetry are non¬ 
singular. The theory is a higher derivative modiflcation of the Einstein action, and 
is obtained by a constructive procedure well motivated in analogy with the analysis 
of point particle motion in special relativity. The resulting theory is asymptotically 
free in a sense which will be specified below. 

Our aim is to construct a theory with the property that the metric approaches 
the de Sitter metric , a metric with maximal symmetry which admits a geodesi- 
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cally complete and nonsingular extension, as the curvature R approaches the Planck 
value Rpi- Here, R stands for any curvature invariant. Naturally, from our classical 
considerations, Rpi is a free parameter. However, if our theory is connected with 
Planck scale physics, we expect Rpi to be set by the Planck scale. 

If successful, the above construction will have some very appealing consequences. 
Consider, for example, a collapsing spatially homogeneous Universe. According to 
Einstein’s theory, this Universe will collapse in hnite proper time to a hnal “big 
crunch” singularity (top left Penrose diagram of Figure 9). In our theory, however, the 
Universe will approach a de Sitter model as the curvature increases. If the Universe is 
closed, there will be a de Sitter bounce followed by re-expansion (bottom left Penrose 
diagram in Figure 9). Similarly, in our theory spherically symmetric vacuum solutions 
would be nonsingular, i.e., black holes would have no singularities in their centers. 
The structure of a large black hole would be unchanged compared to what is predicted 
by Einstein’s theory (top right. Figure 9) outside and even slightly inside the horizon, 
since all curvature invariants are small in those regions. However, for r —0 (where 
r is the radial Schwarzschild coordinate), the solution changes and approaches a de 
Sitter solution (bottom right. Figure 9). This would have interesting consequences 
for the black hole information loss problem. 



FIG. 9. Penrose diagrams for collapsing Universe (left) and black hole (right) in Ein¬ 
stein’s theory (top) and in the nonsingular Universe (bottom). C, E, DS and H stand for 
contracting phase, expanding phase, de Sitter phase and horizon, respectively, and wavy 
lines indicate singularities. 
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To motivate our effective action construction, we turn to a well known analogy, 
point particle motion in the theory of special relativity. 

An Analogy 

The transition from the Newtonian theory of point particle motion to the special 
relativistic theory transforms a theory with no bound on the velocity into one in 
which there is a limiting velocity, the speed of light c Qn the following we use units 
in which fi = c = 1). This transition can be obtained by starting with the action 
of a point particle with world line x{t)\ 

5'oid = J , (44) 

introducing H a Lagrange multiplier held (pwhich couples to the quantity to be 
made hnite, and which has a potential V{^p). The new action is 


Snew / dt 


X + (fx — V{ip) 


(45) 


From the constraint equation 



(46) 


it follows that x^ is limited provided V ((p) increases no faster than linearly in (p for 
large \(p\. The small (p asymptotics of V{(p) is determined by demanding that at low 
velocities the correct Newtonian limit results: 


V{(p) ~ as Iv?! —>• 0 , 
V(p) ^ p as —>• cxo . 


Choosing 


the simple interpolating potential 

2p‘^ 


VW) = 


l + 2f’ 


(47) 

(48) 


(49) 


the Lagrange multiplier can be integrated out, resulting in the well-known action 


Snprw 


1 

2 


/ 


dtVl — x^ 


(50) 


for point particle motion in special relativity. 

Construction 

Our procedure for obtaining a nonsingular Universe theory il is based on gen¬ 
eralizing the above Lagrange multiplier construction to gravity. Starting from the 
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Einstein action, we can introdnce a Lagrange mnltiplier (pi conpled to the Ricci scalar 
R to obtain a theory with limited R\ 

S = J d^x^/^{R + ipiR + Vi{ipi)) , (51) 


where the potential Vi{ipi) satishes the asymptotic conditions (|47|) and (^81). 

However, this action is insnfficient to obtain a nonsingnlar gravity theory. For ex¬ 
ample, singnlar solntions of the Einstein eqnations with i? = 0 are not effected at all. 
The minimal reqnirements for a nonsingnlar theory is that ^ curvatnre invariants 
remain bonnded and the space-time manifold is geodesically complete. Implement¬ 
ing the limiting cnrvatnre hypothesis @, these conditions can be rednced to more 
manageable ones. First, we choose one cnrvatnre invariant and demand that 

it be explicitely bonnded, i.e., |/i| < where if is the Planck scale valne of R. 
In a second step, we demand that as approaches if, the metric g^i, approach 

the de Sitter metric gff, a dehnite nonsingnlar metric with maximal symmetry. In 
this case, all cnrvatnre invariants are antomatically bonnded (they approach their de 
Sitter valnes), and the space-time can be extended to be geodesically complete. 

Onr approach is to implement ^e second step of the above procednre by another 
Lagrange mnltiplier constrnction 1^. We look for a cnrvatnre invariant Rigf^R with 
the property that 

hig,^u) = o ^ = (52) 


introdnce a second Lagrange mnltiplier held (p 2 which conples to R and choose a 
potential V 2 (‘^ 2 ) which forces R to zero at large |(p 2 |: 


S — J Rxy/—g[R + iflR + Vi((pi) -|- (^2-^2 + V2(</52)] , 

with asymptotic conditions (l4^) and (^8]) for VRipi) and conditions 


V 2 {g^ 2 ) ~ const as \ip 2 \ —^ oo 
V2{(P2) ^ aslv^al^o. 


(53) 

(54) 

(55) 


for 1 / 2 ( 952 ). The hrst constraint forces R to zero, the second is reqnired in order to 
obtain the correct low cnrvatnre limit. 

These general conditions are reasonable, bnt not snfhcient in order to obtain a 
nonsingnlar theory. It must still be shown that all solutions are well behaved, i.e., 
that they asymptotically reach the regions |(p 2 | —^ 00 of phase space (or that they 
can be controlled in some other way). This must be done for a specihc realization of 
the above general construction. 

Specific Model 

At the moment we are only able to hnd an invariant R which singles out de 
Sitter space (by demanding R = 0) provided we assume that the metric has special 
symmetries. The choice 
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( 56 ) 




singles out the de Sitter metric among all homogeneous and isotropic metrics (in 
which case adding C^, the Weyl tensor square, is superfluous), all homogeneous and 
anisotropic metrics, and all radially symmetric metrics. 

We choose the action 


S= d^x^J—g 


R + (pii? — ((^2 + ^ 1 ( 9 ^ 1 ) + ^ 2 ( 9 ^ 2 ) 


(57) 


with 




VlHl^ h_Ml±w)) 
“ 1 + v’l V 1 + v’l / 


(58) 


2 

n(^2) = -2^3 4/^-^. (59) 

i + <y92 

The general equations of motion resulting from this action are quite messy. How¬ 
ever, when restricted to homogeneous and isotropic metrics of the form 

ds^ = dd — a{t)‘^{dx^ -|- dy'^ + dz ^), (60) 

the equations are fairly simple. With H = d/a, the two (pi and (^2 constraint equations 
are 

( 61 ) 


H = 


2^/3 


^ 2 , 


(62) 


and the dynamical (^oo equation becomes 


3(1 — ^ 2 ) — {ip2 + 3H(p2) 


(63) 


The phase space of all vacuum configurations is the half plane {(921 > 0, 922 )}- Equa¬ 
tions (|UT|) and (p^ can be used to express H and H in terms of ipi and <^ 22 . The 


(64) 


remaining dynamical equation (|63|) can then be recast as 


d(p2 

dipi 


Vi' 


4V' 


-\i3cp2 + (1 ~ 2 (pi) — i^(hi -|- V2) 


The solutions can be studied analytically in the asymptotic regions and numerically 
throughout the entire phase space. 
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FIG. 10. Phase diagram of the homogeneous and isotropic solutions of the nonsingular 
Universe. The asymptotic regions are labelled by A, B, C and D, flow lines are indicated 
by arrows. 


The resulting phase diagram of vacuum solutions is sketched in Fig. 10 (for 
numerical results, see |^^). The point ( 921 , (p 2 ) = (0,0) corresponds to Minkowski 
space-time M^, the regions I 922 I —^ cxd to de Sitter space. As shown, all solutions 
either are periodic about or else they asymptotically approach de Sitter space. 
Hence, all solutions are nonsingular. This conclusion remains unchanged if we add 
spatial curvature to the model. 

One of the most interesting properties of our theory is asymptotic freedom 0, 
i.e., the coupling between matter and gravity goes to zero at high curvatures. It is 
easy to add matter (e.g., dust or radiation) to our model by taking the combined 
action 


S = Sg + Sm, (65) 

where Sg is the gravity action previously discussed, and Sm is the usual matter action 
in an external background space-time metric. 

We hnd 0 that in the asymptotic de Sitter regions, the trajectories of the solu¬ 
tions in the {(pi, (P 2 ) plane are unchanged by adding matter. This applies, for example, 
in a phase of de Sitter contraction when the matter energy density is increasing expo¬ 
nentially but does not affect the metric. The physical reason for asymptotic freedom 
is obvious: in the asymptotic regions of phase space, the space-time curvature ap¬ 
proaches its maximal value and thus cannot be changed even by adding an arbitrary 
high matter energy density. 

Naturally, the phase space trajectories near ((pi, 922 ) = (0,0) are strongly effected 
by adding matter. In particular, ceases to be a stable fixed point of the evolution 
equations. 
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Discussion 


We have shown that a class of higher derivative extensions of the Einstein theory 
exist for which many interesting solntions are nonsingnlar. Onr class of models is very 
special. Most higher derivative theories of gravity have, in fact, mnch worse singularity 
properties than the Einstein theory. What is special about our class of theories is 
that they are obtained using a well motivated Lagrange multiplier construction which 
implements the limiting curvature hypothesis. We have shown that 

i) all homogeneous and isotropic solutions are nonsingular ii’0 

ii) the two-dimensional black holes are nonsingular B 
hi) nonsingular two-dimensional cosmologies exist 0. 

By construction, all solutions are de Sitter at high curvature. Thus, the theories 
automatically have a ^riod of inflation (driven by the gravity sector in analogy to 
Starobinsky inflation 0) in the early Universe. 

A very important property of our theories is asymptotic freedom. This means that 
the coupling between matter and gravity goes to zero at high curvature, and might 
lead to an automatic suppression mechanism for scalar fluctuations. 


E. Reheating in Inflationary Cosmology 


Reheating is an important stage in inflationary cosmology. It determines the 
state of the Universe after inflation and has consequences for baryogenesis, defect 
formation, and, as will be shown below, maybe even for the composition of the dark 
matter of the Universe. 

After slow rolling, the inflaton held begins to oscillate uniformly in space about 
the true vacuum state. Quantum mechanically, this corresponds to a coherent state 
of /c = 0 inflaton particles. Due to interactions of the inflaton with itself and with 
other helds, the coherent state will decay into quanta of elementary particles. This 
corresponds to post-inflationary particle production. 

Reheating is usually studied using simple scalar held toy models. The one we will 
adopt here consists of two real scalar helds, the inflaton ip with Lagrangian 


]^d^ipd^^ip - ^A((p2 - 


( 66 ) 


interacting with a massless scalar held y representing ordinary matter. The interac¬ 
tion Lagrangian is taken to be 



(67) 


Self interactions of y are neglected. 
By a change of variables 


if = if + a , 
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( 68 ) 


the interaction Lagrangian can be written as 


Cl 


2 ~ 2 I 2 ~2 2 

9 +^9 ^ X ■ 


( 69 ) 


During the phase of coherent oscillations, the held (p oscillates with a frequency 


uj = m^p = 


( 70 ) 


(neglecting the expansion of the Universe which can be taken into account as in 

H,ID- 

Elementary Theory of Reheating 


According to the elementary theory of reheating (see e.g. |^ and |^), the decay 
of the inhaton is calculated using hrst order perturbations theory. According to the 
Feynman rules, the decay rate Vb oi ip (calculated assuming that the cubic coupling 
term dominates) is given by 






( 71 ) 


The decay leads to a decrease in the amplitude of (p (from now on we will drop 
the tilde sign) which can be approximated by adding an extra damping term to the 
equation of motion for ip\ 


(/3 + 3R(p + = -U'((p). (72) 

From the above equation it follows that as long as R > F^, particle production is 
negligible. During the phase of coherent oscillation of (p, the energy density and hence 
H are decreasing. Thus, eventually = F^, and at that point reheating occurs (the 
remaining energy density in (p is very quickly transferred to y particles. 

The temperature Tr at the completion of reheating can be estimated by computing 
the temperature of radiation corresponding to the value of H at which H = Tb- From 
the FRW equations it follows that 

Tr ~ {TBrnpif/^ (73) 

If we now use the “naturalness” constraint]^ 

~ A (74) 


^At one loop order, the cubic interaction term will contribute to A by an amout AA ~ 5^. A 
renormalized value of A smaller than needs to be finely tuned at each order in perturbation 
theory, which is “unnatural”. 
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in conjunction with the constraint on the value of A from (^), it follows that for 

a < iripi, 


Tr < 10^°GeV. 


(75) 


This would imply no GUT baryogenesis, no GUT-scale defect production, and no 
gravitino problems in supersymmetric models with 7 / 13/2 > Tr, where 77 / 3/2 is the 
gravitino mass. As we shall see, these conclusions change radically if we adopt an 
improved analysis of reheating. 


Modern Theory of Reheating 

However, as was hrst realized in |^, the above analysis misses an essential point. 
To see this, we focus on the equation of motion for the matter held y coupled to 
the inhaton (p via the interaction Lagrangian Cj of (^). Taking into account for the 
moment only the cubic interaction term, the equation of motion becomes 


X + ^Hx 


a 


777 , 


= 0. 


( 76 ) 


Since the equation is linear in y, the equations for the Fourier modes Xk decouple: 


Xk + SHxk + (/Cp + 77/J + 2g‘^aip)xk = 0, (77) 

where kp is the time-dependent physical wavenumber. 

Let us for the moment neglect the expansion of the Universe. In this case, the 
friction term in ([77D drops out and kp is time-independent, and Equation ( [77D be¬ 
comes a harmonic oscillator equation with a time-dependent mass determined by the 
dynamics of ip. In the reheating phase, ip is undergoing oscillations. Thus, the mass in 
([TtI) is varying periodically. In the mathematics literature, this equation is called the 
Mathieu equation. It is well known that there is an instability. In physics, the effect 
is known as parametric resonance (see e.g. [^). At frequencies a;„ corresponding 
to half integer multiples of the frequency u of the variation of the mass, i.e. 

T) 

^l = kl + ml = {-u:f 77 = 1,2,..., (78) 

there are instability bands with widths Aa;,i. For values of ujk within the instability 
band, the value of Xk increases exponentially: 

2 

Xk ~ with p ~ ^ , (79) 

CO 

with ipo being the amplitude of the oscillation of ip. Since the widths of the instability 
bands decrease as a power of the (small) coupling constant with increasing n, for 
practical purposes only the lowest instability band is important. Its width is 

Auk ~ • (80) 
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Note, in particular, that there is no ultraviolet divergence in computing the total 
energy transfer from the ip to the x S^ld due to parametric resonance. 

It is easy to include the effects of the expansion of the Universe (see e.g. P8| , |i^Ji5[| ). 
The main effect is that the value of ujk becomes time-dependent. Thus, a mode 
slowly enters and leaves the resonance bands. As a consequence, any mode lies in 
the resonance band for only a hnite time. This implies that the calculation of energy 
transfer is perfectly well-behaved. No inhnite time divergences arise. 

It is now possible to estimate the rate of energy transfer, whose order of magnitude 
is given by the phase space volume of the lowest instability band multiplied by the 
rate of growth of the mode function Xk- Using as an initial condition for Xk the value 
Xk ^ H given by the magnitude of the expected quantum fluctuations, we obtain 




From (^TD it follows that provided that the condition 


pAt >> 1 


(81) 


(82) 


is satished, where At < H~^ is the time a mode spends in the instability band, 
then the energy transfer will precede fast on the time scale of the expansion of the 
Universe. In this case, there will be explosive particle production, and the energy 
density in matter at the end of reheating will be given by the energy density at the 
end of inflation. 

The above is a summary of the main physics of the modern theory of reheating. 
The actual analysis can be rehned in many ways (see e.g. |]4^ , ^ ). First of all, it 
is easy to take the expansion of the Universe into account explicitly (by means of a 
transformation of variables), to employ an exact solution of the background model 
and to reduce the mode equation for Xk to a Hill equation, an equation similar to the 
Mathieu equation which also admits exponential instabilities. 

The next improvement consists of treating the x held quantum mechanically (keep¬ 
ing ip as & classical background held). At this point, the techniques of quantum held 
theory in a curved background can be applied. There is no need to impose artihcial 
classical initial conditions for Xk- Instead, we may assume that y starts in its initial 


vacuum state (excitation of an initial thermal state has been studied in ||^), and the 


Bogoliubov mode mixing technique (see e.g. |^) can be used to compute the number 
of particles at late times. 

Using this improved analysis, we recover the result dSlf) . Thus, provided that the 
condition (p^ is satished, reheating will be explosive. Working out the time At that 
a mode remains in the instability band for our model, expressing H in terms of ipQ 
and rripi, and uj in terms of a, and using the naturalness relation X, the condition 
for explosive particle production becomes 


ipoTflpi 


>> 1 


(83) 





















which is satished for all chaotic inflation models with a < nipi (recall that slow rolling 
ends when ip ~ mpi and that therefore the initial amplitude (po of oscillation is of the 
order rripi). 

We conclude that rather generically, reheating in chaotic inflation models will be 
explosive. This implies that the energy density after reheating will be approximately 
equal to the energy density at the end of the slow rolling period. Therefore, as 
suggested in and [^, respectively, GUT scale defects may be produced after 

reheating and GUT-scale baryogenesis scenarios may be realized, provided that the 
GUT energy scale is lower than the energy scale at the end of slow rolling. 

Note, however, that the state of x after parametric resonance is not a thermal 
state. The spectrum consists of high peaks in distinct wave bands. An important 
question which remains to be studied is how this state thermalizes. For some inter¬ 
esting work on this issue see |^. As emphasized in and [Q, the large peaks 


in the spectrum may lead to symmetry restoration and to the efficient production of 
topological defects (for a differing view on this issue see [^,^]). Since the state after 
explosive particle production is not a thermal state, it is useful to follow and call 
this process “preheating” instead of reheating. 

A father interesting conjecture which emerges from the parametric resonance anal¬ 
ysis of preheating cSB is that the dark matter in the Universe may consist of remnant 
coherent oscillations of the inflaton field. In fact, it can easily be checked from (^^ 
that the condition for efficient transfer of energy eventually breaks down when (po 
has decreased to a sufficiently small value. For the model considered here, an order 
of magnitude calculation shows that the remnant oscillations may well contribute 
significantly to the present value of G. 

Note that the details of the analysis of preheating are quite model-dependent. 
In fact 0, in many models one does not get the kind of “narrow-band” resonance 
discussed here, but “wide-band” resonance. In this case, the energy transfer is even 
more efficient. 

There has recently been a lot of work on various aspects of reheating (see e.g. 
|]58H^ for different approaches). Many important questions, e.g. concerning ther- 
malization and back-reaction effects during and after preheating (or parametric res¬ 
onance) remain to be fully analyzed. 


F. Summary 


The inflationary Universe is an attractive scenario for early Universe cosmology. 
It can resolve some of the problems of standard cosmology, and in addition gives rise 
to a predictive theory of structure formation (see e.g. [^| for a recent review). 

However, important unsolved problems of principle remain. Rather generically, 
the predicted amplitude of perturbations is too large (the spectral shape, however, is 
in quite good agreement with the observations). The present realizations of inflation 
based on scalar field also make the cosmological constant problem more accute. In 
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addition, there are no convincing particle-physics based realizations of inflation. Many 
models of inflation resort to introducing a new matter sector. It is important to search 
for a better connection between modern particle physics / held theory and inhation. 
String cosmology and dilaton gravity (see e.g. the recent reviews in [Q) may provide 
an interesting new approach to the unihcation of inhation and fundamental physics. 

Recently, there has been much progress in the understanding of the energy transfer 
at the end of inhation between the inhaton held and matter. It appears that resonance 
phenomena such as parametric resonance play a crucial role. These new reheating 
scenarios lead to a high reheating temperature, although much more work remains to 
be done before one can reach a hnal conclusion on this issue. 


III. CLASSICAL AND QUANTUM THEORY OF COSMOLOGICAL 

PERTURBATIONS 

In inhationary Universe and topological defect models of structure formation, 
small amplitude seed perturbations are predicted to arise due to particle physics 
ehects in the very early Universe. They then grow by gravitational instability to 
produce the cosmological structures we observe today. In order to be able to make 
the connection between particle physics and observations, it is important to under¬ 
stand the gravitational evolution of huctuations. This section will introduce the basic 
concepts of this topic. 

As is evident from Figure 5 and from the discussion of inhation in the previous 
section, general relativity and quantum mechanics both play a fundamental role in 
the theory of perturbations. In inhationary Universe models, quantum ehects seed 
the huctuations, and thus a quantum analysis of the generation of huctuations is 
essential. However, since the huctuations are small, a linearized analysis is sufficient. 
Since the scales on which we are interested in following the huctuations are larger than 
the Hubble radius for a long time interval, Newtonian gravity is obviously inadequate 
to treat these perturbations, and general relativistic ehects become essential. 

In this section, we will hrst introduce some basic notation, then discuss the New¬ 
tonian theory of linear huctuations before turning to the full relativistic analysis. 


A. Basic Issues 

In this article we only discuss theories in which structures grow by gravitational 
accretion. The basic mechanism is easy to understand. Consider hrst a hat space- 
time background. A density perturbation with 6p > 0 leads to an excess gravitational 
attractive force F acting on the surrounding matter. This force is proportional to 6p, 
and will hence lead to exponential growth of the perturbation since 

6p rsj P rsj 6p ^ 6p exp(Q;t) (84) 
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with a constant a which is proportional to Newton’s constant G. 

In an expanding background space-time, the acceleration is damped by the ex¬ 
pansion. If r{t) is the physical distance of a test particle from the perturbation, then 
on a scale r 


5p 


6p 

r2(t) ’ 


( 86 ) 


which results in power-law increase of Sp. The goal of this subsection is to discuss the 
growth rates of inhomogeneities in more detail (see e.g. |6^ , |65[| for modern reviews). 


Because of our assumption that all perturbations start out with a small amplitude, 
we can linearize the equations for gravitational fluctuations. The analysis is then 
greatly simplihed by going to momentum space in which all modes 6{k) decouple. 
We expand the fractional density contrast 6{x) as follows: 


= (27r)-3/Vi/2 f d^ke^^-^6{k) , 
p J 

where Id is a cutoff volume which disappears from all physical observables. 
The “power spectrum” P{k) is dehned by 


( 86 ) 


P(k) =< |6(*:)P > 


(87) 


where the braces denote an ensemble average (in most structure formation models, 
the generation of perturbations is a stochastic process, and hence observables can 
only be calculated by averaging over the ensemble. For observations, the braces can 
be viewed as an angular average). 

The physical measure of mass fluctuations on a length scale A is the r.m.s. mass 
fluctuation 6M/M{X) on this scale. It is determined by the power spectrum in the 
following way. We pick a center Xg of a sphere Bx{x^) of radius A and calculate 


5M 

W 





d^x6(x) 


1 |2 

V(B ,)' 


( 88 ) 


where V{B\) is the volume of the sphere. Inserting the Fourier decomposition (^^ 
and taking the average value of this quantity over all Xg yields 


< 


'6M" 


(A) >= /dW,(A)|5(fc)p 


with a window function Wk{X) with the following properties 

Wk{X) 

Therefore the r.m.s. mass perturbation on a scale A becomes 


~ 1 k < kx = 2 tt/X 
~ 0 k > kx ■ 


(89) 


(90) 
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< >-klP{k,). (91) 

If P{k) ~ fc" then n is called the index of the power spectrum. For n = 1 we get 
the so-called Harrison-Zel’dovich scale invariant spectrum 

Both inflationary Universe and topological defect models of structure formation 
predict a roughly scale invariant spectrum. The distinguishing feature of this spec¬ 
trum is that the r.m.s. mass perturbations are independent of the scale k when 
measured at the time tnik) when the associated wavelength is equal to the Hubble 
radius, i.e., when the scale “enters” the Hubble radius. Let us derive this fact for the 
scales entering during the matter dominated epoch. The time tnik) is determined by 

k~^a{tHik)) = tnik) (92) 


which leads to tnik) ~ k~^. According to the linear theory of cosmological perturba¬ 
tions discussed in the following subsection, the mass fluctuations increase as a(t) for 
t > teq- Hence 


6M 


{k,tH{k)) 




6M 


{k,t) const, 


(93) 


since the hrst factor scales (from (p^ as k and - using (^) and inserting n = 1 - 
the second as 


B. Newtonian Theory 

The Newtonian theory of cosmological perturbations is an approximate analysis 
which is valid on wavelengths A much smaller than the Hubble radius t and for 
negligible pressure p, i.e., p ^ p. It is based on expanding the hydrodynamical 
equations about a homogeneous background solution. 

The starting points are the continuity, Euler and Poisson equations 

P + V(pi;) = 0 (94) 


L + (t ■ V)i; = — V0-Vp 

P 


(96) 


VV = 4!rGp (96) 

for a fluid with energy density p, pressure p, velocity v and Newtonian gravitational 
potential 0, written in terms of physical coordinates {t, r). 

The transition to an expanding space is made by introducing comoving coordinates 
X and peculiar velocity u = P. 
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(97) 


r = a(t)x 


v = a{t)x + a{t)u. (98) 

The first term on the right hand side of (|^) is the expansion velocity. 

The perturbation equations are obtained by linearizing Equations - 0) about 
a homogeneous background solution p = pit), p = 0 and u = 0. The linearization 
ansatz can be written 

p{x,t) = p{t){l + 6{x,t)) 1^1 << 1. (99) 

If we consider adiabatic perturbations (no entropy density variations), then after some 
algebra the linearized equations become 

5 + V-m = 0, (100) 

u + 2\ =-a^(V6(j) + clV6) ( 101 ) 

a 

and 


= 47rGpa^6 , 

with the speed of sound Cg given by 




dp 

dp 


( 102 ) 


(103) 


The two first order equations d® and m can be combined to yield a single second 
order differential equation for 6. With the help of ( p.02|) this equation reads 


5 + 2H5 - AirGpd 

which in momentum space becomes 


v72 


= 0 




6k + ‘2H5k + 


AirGp ) = 0 


(104) 


(105) 


Here, H{t) as usual denotes the expansion rate, and 5k stands for 5{k). 

Already a quick look at Equation (|105|) reveals the presence of a distinguished 


(106) 


scale for cosmological perturbations, the Jeans length 

271 


A.7 = 


kj 
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with 


2 



AirGp 


(107) 


On length scales larger than Aj, the spatial gradient term is negligible, and the term 
linear in S in (11051) acts like a negative mass square quadratic potential with damping 
due to the expansion of the Universe, in agreement with the intuitive analysis leading 
to (refintul) and (^). On length scales smaller than Aj, however, (|105|) becomes a 
damped harmonic oscillator equation and perturbations on these scales decay. 

For t > tf>g and for A Aj, Equation (|105|) becomes 

h + = 0 ( 108 ) 

and has the general solution 

Sk{t) = + C2t~^. (109) 


This demonstrates that for t > ten and A S> Aj, the dominant mode of perturbations 


increases as a(f), a result we already used in the previous subsection (see 
For A <C Aj and t > teq. Equation (|105|) becomes 


!))• 


'Sk + “iHdk + ~ 0 , 

and has solutions corresponding to damped oscillations: 

Sk{t) ~ a“^'^^(f) exp{±ics/c J dt'a{t')~^} . 


( 110 ) 


( 111 ) 


As an important application of the Newtonian theory of cosmological perturba¬ 
tions, let us compare sub-horizon scale fluctuations in a baryon-dominated Universe 
(fl = = 1) and in a CDM-dominated Universe with flcDM = 0.9 and hi = 1. We 

consider scales which enter the Hubble radius at about teg. 

In the initial time interval teg < t < tree, the baryons are coupled to the photons. 
Hence, the baryonic fluid has a large pressure pb 

1 , , 

Pb ^ Pr = -Pr , (112) 

and therefore the speed of sound is relativistic 


Cs 





1 

73 



(113) 
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The value of Cs slowly decreases in this time interval, attaining a value of about 1/10 
at tree- From ( 107 ) it follows that the Jeans mass Mj, the mass inside a sphere of 


radius Aj, increases until tree when it reaches its maximal value 


A-TT 

Mr" = Mj((,„) = — ~ 


(114) 


At the time of recombination, the baryons decouple from the radiation fluid. 
Hence, the baryon pressure pb drops abruptly, as does the Jeans length (see ( |107|) ). 
The remaining pressure ps is determined by the temperature and thus continues to 
decrease as t increases. It can be shown that the Jeans mass continues to decrease 
after tree, starting from a value 




(115) 


(where the superscript ” indicates the mass immediately after teq. 

In contrast, CDM has negligible pressure throughout the period t > teq and hence 
experiences no Jeans damping. A CDM perturbation which enters the Hubble radius 
at teq with amplitude 5i has an amplitude at tree given by 

4°“'((„,)c. ^^4.. (116) 

whereas a perturbation with the same scale and initial amplitude in a baryon- 
dominated Universe is damped 

(^i ‘'A.. (117) 

In order for the perturbations to have the same amplitude today, the initial size of 
the inhomogeneity must be much larger in a BDM-dominated Universe than in a 
CDM-dominated one: 

4fo«((,,) M ■ (118) 

For D = 1 and h = l/2 the enhancement factor is about 30. 

In a CDM-dominated Universe the baryons experience Jeans damping, but after 
tree the baryous quickly fall into the potential wells created by the CDM perturbations, 
and hence the baryon perturbations are proportional to the CDM inhomogeneities. 

The above considerations, coupled with information about CMB anisotropies, can 
be used to rule out a model with D = = 1. The argument goes as follows. For 

adiabatic fluctuations, the amplitude of CMB anisotropies on an angular scale d is 
determined by the value of 5p/p (strictly speaking, the relativistic potential <h to be 
discussed in the following subsection) on the corresponding length scale A('d) at teq'. 
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(119) 


On scales of clusters we know that (for 0 = 1 and h = 1/2) 

f—) teg) ^ z{teg)~^ ~ 10"^ , (120) 

\ P / CDM 


using the fact that today on cluster scales 5p/p ~ 1. The bounds on 6T/T on small 
angular scales are 


~T 


( d ) <<10 ^, 


( 121 ) 


consistent with the predictions for a CDM model, but inconsistent with those of a 
O = Qb = 1 model, according to which we would expect anisotropies of the order of 
10“^. This is yet another argument in support of the existence of nonbaryonic dark 
matter. 

To conclude this subsection, let us briefly discuss two further aspects related 
to Newtonian perturbations. The first concerns matter inhomogeneities during the 
radiation-dominated epoch. We consider matter fluctuations with = 0 in a smooth 
relativistic background. In this case. Equation (|105|) becomes 


Sk + ^HSk — AirCpmSk — 0 , 


( 122 ) 


where pm denotes the average matter energy density. The Hubble expansion param¬ 
eter obeys 




SttG 


(Pm T Pr) 5 


(123) 


with pr the background radiation energy density. For t < teq, Pm is negligible in both 
( p.22|) and (|123|) , and (|122|) has the general solution 


St(t) = Cl log! + C2. 


(124) 


In particular, this result implies that CDM perturbations which enter the Hubble 
radius before teg have an amplitude which grows only logarithmically in time until 

teg- 

Finally, we consider hot dark matter (HDM) fluctuations. Whereas CDM particles 
are cold, i.e. their peculiar velocity is negligible for all times relevant for structure 
formation, HDM particles have relativistic velocities at teg, i.e. v{teq) ~ 1. The prime 
candidate for HDM is a 25/i5QeV tan neutrino. 

The new aspect of HDM is related to neutrino free streaming 0. Because of 
the large velocity of the dark matter particles, pure dark matter inhomogeneities are 
washed out on all scales below the neutrino free streaming length Xj{t), 
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(126) 


A‘(«) ~ v{t)z{t)t, 

which is the comoving distance the particles move in one Hubble expansion time. 
Since the neutrino velocity v(t) and the redshift z(t) both scale as the free 

streaming length decreases as 

~ (126) 

after teq (before teg the radiation pressure dominates). 

Hence, in an inflationary HDM model in which the fluctuations are dark matter 
inhomogeneities, all perturbations on scales A smaller than the maximal value of A^(f) 
are erased. The critical scale is given by the value of A^(f) at the time when the 
neutrinos become non-relativistic, which is in turn determined by the neutrino mass 
ruy. The result is 

A»“^30Mpc(^)"A (127) 

a scale much larger than the mean separation of galaxies and clusters. Since we ob¬ 
serve galaxies outside of large-scale structures, this model is in blatant disagreement 
with observations. However, theories in which the primordial perturbations are nona- 
diabatic long-lived seeds (e.g. cosmic strings), may well be viable if the dark matter 
is hot. As we shall see in Section 4, the cosmic string model in fact works well for 
hot dark matter. 


C. Relativistic Theory: Classical Analysis 

On scales larger than the Hubble radius (A > t) the Newtonian theory of cos¬ 
mological perturbations obviously is inapplicable, and a general relativistic analysis 
is needed. On these scales, matter is essentially frozen in comoving coordinates. 
However, space-time fluctuations can still increase in amplitude. 

In principle, it ^ straightforward to work out the general relativistic theory of 
linear fluctuations S. We linearize the Einstein equations 

= 8nGT^, (128) 

(where Gq_y is the Einstein tensor associated with the space-time metric and 
is the energy-momentum tensor of matter) about an expanding FRW background 

t) = + h^^{x, t) (129) 

t) = t) (130) 

and pick out the terms linear in hq,y and 5(p to obtain 
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SGf,, = MG5T^,. (131) 

In the above, is the perturbation in the metric and is the fluctuation of the 
matter held ip. We have denoted all matter fields collectively by ip. 

In practice, there are many complications which make this analysis highly nontriv¬ 
ial. The first problem is “gauge invariance” ® Imagine starting with a homogeneous 
FRW cosmology and introducing new coordinates which mix a; and t. In terms of the 
new coordinates, the metric now looks inhomogeneous. The inhomogeneous piece of 
the metric, however, must be a pure coordinate (or ’’gauge”) artefact. Thus, when 
analyzing relativistic perturbations, care must be taken to factor out effects due to 
coordinate transformations. 



FIG. 11. Sketch of how two choices of the mapping from the background space-time 
manifold Alo to the physical manifold A1 induce two different coordinate systems on Al. 


The issue of gauge dependence is illustrated in Fig. 11. A coordinate system on 
the physical inhomogeneous space-time manifold Ad can be viewed as a mapping T> 
of an unperturbed space-time Ado into Ad. A physical quantity Q is a geometrical 
function defined on Ad. There is a corresponding physical quantity defined on 
Ado. In the coordinate system given by V, the perturbation 5Q of Q at the space-time 
point p e Ad is 


SQ(.P)=Q{p)-'°P{D-'(p)). (132) 

However, in a second coordinate system V the perturbation is given by 

SQ{p)^Q{p)-i"'Q{V-'{p)). (133) 


The difference 
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AQ{p) = SQ(p) - iQ(p) 


(134) 


is obviously a gauge artefact and carries no physical meaning. 

There are various methods of dealing with gauge artefacts. The simplest and 
most physical approach is to focus on gauge invariant variables, i.e., combinations 
of the metric and matter perturbations which are invariant under linear coordinate 
transformations. 

The gauge invariant theory of cosmological perturbations is in principle straight¬ 
forward, although technically rather tedious. In the following I will summarize the 
main steps and refer the reader to [^] for the details and further references (see also 
0 for a pedagogical introduction and ||7T|- [7^] for other approaches). 

We consider perturbations about a spatially flat Friedmann-Robertson-Walker 
metric 


ds^ = a^{ri){drf‘ — dx^) 


(135) 


where rj is conformal time (related to cosmic time t by a{ri)dr] = dt). A scalar metric 
perturbation (see for a precise dehnition) can be written in terms of four free 


functions of space and time: 


= a^iv) 



\ 

2('ipSij + E,ij)) 


(136) 


Scalar metric perturbations are the only perturbations which couple to energy density 
and pressure. 

The next step is to consider inhnitesimal coordinate transformations 


x^^' = 


(137) 


which preserve the scalar nature of and to calculate the induced transformations 
of and E. Then we hnd invariant combinations to linear order. (Note that 

there are in general no combinations which are invariant to all orders S.) After some 
algebra, it follows that 


$ = (j) + a-^[{B 


T = V’ 


(B 


- E')a]' 
E') 


(138) 

(139) 


are two invariant combinations. In the above, a prime denotes differentiation with 
respect to rj. 

There are various methods to derive the equations of motion for gauge invariant 
variables. Perhaps the simplest way 0 is to consider the linearized Einstein equations 


( p.31|) and to write them out in the longitudinal gauge dehned by 


B = E = 0 


(140) 
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and in which $ = 0 and d' = 0, to directly obtain gange invariant eqnations. 

For several types of matter, in particnlar for scalar held matter, the pertnrbation 
of has the special property 




(141) 


which imples <F = d'. Hence, the scalar-type cosmological pertnrbations can in this 
case be d^ri^^by a single gange invariant variable. The eqnation of motion takes 
the form 




(142) 


where 


£ =-h <h 

^ 3 1 + w; 


(143) 


The variable w = p/p (with p and p backgronnd pressnre and energy density 
respectively) is a measnre of the backgronnd eqnation of state. In particnlar, on 
scales larger than the Hnbble radins, the right hand side of (|142|) is negligible, and 
hence £ is constant. 

The resnlt that £ = 0 is a very powerfnl one. Let ns hrst imagine that the eqnation 
of state of matter is constant, he., w = const. In this case, £ = 0 implies 


<h(t) = const, 


(144) 


he., this gange invariant measnre of pertnrbations remains constant ontside the Hnb¬ 
ble radins. 

Next, consider the evolntion of <F dnring a phase transition from an initial phase 
with w = Wi io a. phase with w = Wf. Long before and after the transition, <F is 
constant becanse of (|144|) , and hence £ = 0 becomes 




1 -I- te 


<h = const, 


(146) 


In order to make contact with matter pertnrbations and Newtonian intnition, it 
is important to remark that, as a conseqnence of the Einstein constraint eqnations, 
at Hnbble radius crossing $ is a measure of the fractional density fluctuations: 

<F(/c, tnik)) ~ — {k, tnik)). (146) 

P 

(Note that the latter quantity is approximately gauge invariant on scales smaller than 
the Hubble radius). 
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D. Relativistic Theory: Quantum Analysis 


The question of the origin of classical density perturbations from quantum fluc¬ 
tuations in the de Sitter phase of an inflationary Universe is a rather subtle issue. 
Starting from a homogeneous quantum state {e.g., the vacuum state in the FRW 
coordinate frame at time tj, the beginning of inflation), a naive semiclassical anaylsis 
would predict the absence of fluctuations since < > is independent of 

space. 

However, as a simple thought experiment shows, such a naive analysis is inap¬ 
propriate. Imagine a local gravitational mass detector D positioned close to a large 
mass M which is suspended from a pole. The decay of an alpha particle will sever 
the cord by which the mass is held to the pole and the mass will drop. According to 
the semiclassical prescription 

= SttG < , (147) 

the metric {i.e., the mass measured) will slowly decrease. This is obviously not what 
happens. The mass detector shows a signal which corresponds to one of the classical 
trajectories which make up the state \'ip >, a trajectory corresponding to a sudden 
drop in the gravitational force measured. 

The origin of classical density perturbations as a consequence of quantum fluctua¬ 
tions in a homogeneous state \ip > can be analyzed along similar lines. The quantum 
to classical transition is picking out one of the typical classical trajectories 

which make up the wave function of \'ijj >. We can implement @0 the procedure as 
follows: Define a classical scalar field 

t) = ipo{t) + 6ip{x, t) (148) 

with vanishing spatial average of 6(p. The induced classical energy momentum tensor 
T^l {x, t) which is the source for the metric is given by 

t) = T^^((^d(A, t )), (149) 

where {g>d{x,t)) is defined as the canonical energy-momentum tensor of the clas¬ 
sical scalar field (pci{x,t). Unless 6ip vanishes, T^l is inhomogeneous. 

For applications to chaotic inflation, we take > to be a Gaussian state with 
mean value (po{t) 

<AA‘^{x,t)\'ijj>=ipl{t). (150) 

Its width is taken to be the width of the vacuum state of the free scalar field theory 
with mass determined by the curvature of at (po- This state is used to define 

the Fourier transform 6(p{k,t) by 

\5(p{k)\^ =< %l)\ \Cp{k)\^\A > • (151) 
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The amplitude of 6(p{k) is identified with the width of the ground state wave function 
of the harmonic osciallator (p{k). (Recall that each Fourier mode of a free scalar held 
is a harmonic oscillator). Note that no divergences arise in the above construction. 

By linearizing (|149|) about (^o(^) we obtain the perturbation of the energy- 
momentum tensor. In particular, the energy density huctuation 5p{k) is given by 

6p{k) = pQ5(p{k) + V'{’pQ)6(p{k). (152) 


To obtain the initial amplitude (^9|) of 6M/M, the above is to be evaluated at the 
time ti{k). 

The computation of the spectrum of density perturbations produced in the de 
Sitter phase has been reduced to the evaluation of the expectation value (|151|). First, 


we must specify the state IV' >• (Recall that in non-Minkow^i space-times there is 
no uniquely dehned vacuum state of a quantum held theory Eli). We pick the FRW 
frame of the pre-inhationary period. In this frame, the number density of particles 
decreases exponentially. Hence we choose \'ip > to be the vacuum state in this frame 
for a discussion of other choices). '0[<^(fc),t], the wave functional of \^lJ >, 


see 


can be calculated explicitly. It is basically the superposition of the ground state wave 
functions for all oscillators 


^] =-^exp | —^(27r) j d^koj{k^t)\^p{^[ 

N is a. normalization constant and uj{k,t) ~ FT at t = ti{k). Hence 

6(p{k,t) = (27r)^/^a“^/^a;(fc, ~ (27r)^/^fc“^/^iF, t = ti{k) . 


(153) 


(154) 


Given the above determination of the intitial amplitude of density perturbations 
at the time when they leave the Hubble radius during the de Sitter phase, and the 
general relativistic analysis of the evolution of huctuations discussed in the previous 
subsection, it is easy to evaluate the r.m.s. inhomogeneities when they reenter the 
Hubble radius at time t/(fc). 

First, we combine (^, @, (^), (0) and m to obtain 


If the background scalar field is rolling slowly, then 


'V'MH' 


Po 


and 


(1 + p/p)(ti(k)) ~ . 


(155) 


(156) 


(157) 


Combining (155), 


, (p.57|) and (^ we get 
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5M 


{k, tf{k)) 




3H^ 


^o{U{k)) |0o(ti(A;))| 


(158) 


This result can now be evaluated for specific models of inflation to find the conditions 
on the particle physics parameters which give a value 


6M 

~W 


{k, tf{k)) ~ 10 


-5 


(159) 


which is required if quantum fluctuations from inflation are to provide the seeds for 
galaxy formation and agree with the CMB anisotropy limits. 

For chaotic inflation with a potential 


V((p) = 

we can solve the slow rolling equations for the inflaton to obtain 

6M ,,,, m 

—(kMk)) ~10— 


(160) 


(161) 


which implies that m 10^^ GeV to agree with (|159|) 
Similarly, for a quartic potential 


we obtain 


6M 

~W 


V{^) = -A(p^ 


(fc, tf{k)) ~ 10 ■ 


(162) 


(163) 


which requires A < 10 in order not to conflict with observations. 


The conditions (|161| ) and (|163| ) require the presence of small parameters in the 
particle physics model. It has been shown quite generally 0 that small parameters 
are required if inflation is to solve the fluctuation problem. 

I have chosen to present the analysis of fluctuations in inflationary cosmology in 
two separate steps in order to highlight the crucial physics issues. Having done this, 
it is possible to step back and construct a unihed analysis of the quantum generation 
and classical evolution of perturbations in an inflationary Universe (for a detailed 
review see 
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The basic point is that at the linearized level, the equations describing both grav¬ 
itational and matter perturbations can be quantized in a consistent way. The use of 
gauge invariant variables makes the analysis both physically clear and computation¬ 
ally simple. 

The hrst step of this analysis is to consider the action for the linear perturbations 
in a background homogeneous and isotropic Universe, and to express the result in 
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terms of gauge invariant variables describing the fluctuations. Focusing on the scalar 
perturbations, it turns out that after a lot of algebra the action reduces to the action 
of a single gauge invariant free scalar held with a time dependent mass 110 (the time 
dependence relects the expansion of the background space-time). This result is not 
suprising. Based on the study of classical cosmological perturbations, we know that 
there is only one held degree of freedom for the scalar perturbations (for matter the¬ 
ories which obey (|141[ )). Since at the linearized level there are no mode interactions, 
the action for this held must be that of a free scalar held. 

The action thus has the same form as the action for a scalar matter held in a time 
dependent gravitational or electromagnetic background, and we can use standard 
methods to quantize this theory (see e.g. |]5^). If we employ canonical quantization, 
then the mode functions of the held operator obey the same classical equations as we 
derived in the gauge-invariant analysis of relativistic perturbations. 

The time dependence of the mass is relected in the nontrivial form of the solutions 
of the mode equations. The mode equations have growing modes which correspond 
to particle production or equivalently to the generation and amplihcation of huctua- 
tions. We can start the system oh (e.g. at the beginning of inhation) in the vacuum 
state (dehned as a state with no particles with respect to a local comoving observer). 
The state dehned this way will not be the vacuum state from the point of view of 
an observer at a later time. The Bogoliubov mode mixing technique can be used to 
calculate the number density of perticles at a later time. In particular, expectation 
values of held operators such as the power spectrum can be computed, and the re¬ 
sults agree with those obtained using the heuristic analysis at the beginning of this 
subsection. 


E. Summary 

To summarize the main results of the analysis of density huctuations in inhationary 
cosmology: 

1. Quantum vacuum huctuations in the de Sitter phase of an inhationary Universe 
are the source of perturbations. 

2. The quantum perturbations decohere on scales outside the Hubble radius and 
can hence be treated classically. 

3. As a consequence of the change in the background equation of state, the clas¬ 
sical evolution outside the Hubble radius produces a large amplihcation of the 
perturbations. In fact, unless the particle physics model contains very small 
coupling constants, the predicted huctuations are in excess of those allowed by 
the bounds on cosmic microwave anisotropies. 
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4. The quantum generation and classical evolution of fluctuations can be treated 
in a unified manner. The formalism is no more complicated that the study of a 
free scalar field in a time dependent background. 

5. Inflationary Universe models generically produce an approximately scale invari¬ 
ant Harrison-Zel’dovich spectrum 

~ const. (164) 


It is not hard to construct models which give a different spectrum. All that is 
required is a signihcant change in H during the period of inflation. 

Once inside the Hubble radius, the evolution of the mass perturbations is influ¬ 
enced by the damping effects discussed in Section 3.B, which in turn depend on the 
composition of the dark matter. 

On scales which enter the Hubble radius before t^g, the perturbations can only 
grow logarithmically in time between tf{k) and teg. This implies that (up to loga¬ 
rithmic corrections), the mass perturbation spectrum is flat for wavelengths smaller 
than \eg, the comoving Hubble radius at tf>g: 

^(A,f) ~ const, t<teq,X<Xeq, (165) 

whereas on larger scales 


6M 


(A, t) oc A 


(166) 


Equations (|165|) and (|166|) give the power spectrum in an O = 1 inflationary CDM 
model. If the dark matter is hot, then neutrino free streaming cuts off the power 
spectrum at A™“* (see (|127|) ). 


IV. TOPOLOGICAL DEFECTS, STRUCTURE FORMATION AND 

BARYOGENESIS 

A. Quantifying Data on Large-Scale Structnre 

It is length scales corresponding to galaxies and larger which are of greatest inter¬ 
est in cosmology when attempting to find an imprint of the primordial fluctuations 
produced by particle physics. On these scales, gravitational effects are assumed to 
be dominant, and the fluctuations are not too far from the linear regime. On smaller 
scales, nonlinear gravitational and hydrodynamical effects determine the hnal state 
and mask the initial perturbations. 
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To set the scales, consider the mean separation of galaxies, which is abont 
Mpc 0, and that of Abell clnsters which is aronnd 25h“^ Mpc 0. The largest 
coherent strnctnres seen in cnrrent redshift snrveys have a length of abont 100 h“^ 
Mpc i’i , the recent detections of CMB anisotropies probe the density held on length 
scales of abont Mpc, and the present horizon corresponds to a distance of 

abont 3 ■ Mpc. 

Galaxies are gravitationally bonnd systems containing billions of stars. They are 
non-randomly distribnted in space. A qnantitative measnre of this non-randomness 
is the “two-point correlation fnnction” .^ 2 (^) which gives the excess probability of 
hnding a galaxy at a distance r from a given galaxy: 


=< 


n{r) — no 
no 


> . 


(167) 


Here, no is the average nnmber density of galaxies, and n{r) is the density of galaxies 
a distance r from a given one. The pointed braces stand for ensemble averaging. 

Recent observational resnlts from a various galaxy redshift surveys yield reason¬ 
able agreement @ with a form 

^ (168) 


with scaling length tq — 5h~^ Mpc and power 7 ~ 1.8. A theory of structure 
formation must explain both the amplitude and the slope of this correlation function. 

On scales larger than galaxies, a better way to quantify structure is by means of 
large-scale systematic redshift surveys. Such surveys have discovered coherent planar 
structures and voids on scales of up to 100h“^ Mpc. Fig. 12 is a sketch of redshift 
versus angle a in the sky for one 6 ° slice of the sky 0. The second direction in 
the sky has been projected onto the a — z plane. The most prominent feature is 
the band of galaxies at a distance of about 100h“^ Mpc. This band also appears in 
neighboring slices and is therefore presumably part of a planar density enhancement 
of comoving planar size of at least (50 x 100) x h~‘^ Mpc^. This structure is often 
called the “Great Wall.” The challenge for theories of structure formation is not only 
to explain the fact that galaxies are nonrandomly distributed, but also to predict 
both the observed scale and topology of the galaxy distribution. Topological defect 
models of structure formation attempt to address these questions. 

Until 1992 there was little evidence for any convergence of the galaxy distribution 
towards homogeneity. Each new survey led to the discovery of new coherent structures 
in the Universe on a scale comparable to that of the survey. In 1996, results of a much 
deeper redshift survey were published fl which for the hrst time hnd no coherent 
structures on the scale of the entire survey. In fact, so coherent structures on scales 
larger than 100h“^ Mpc are seen. This is the hrst direct evidence for the cosmological 
principle from optical surveys (the isotropy of the GMB has for a long time been a 
strong point in its support). 
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Right ascension 



FIG. 12. Results from the CFA redshift survey. Radial distance gives the redshift of 
galaxies, the angular distance corresponds to right ascension. The results from several slices 
of the sky (at different declinations) have been projected into the same cone. 


In summary, a lot of data from optical and infrared galaxies are currently available, 
and new data are being collected at a rapid rate. The observational constraints on 
theories of structure formation are becoming tighter. 


B. Toplogical Defects 

According to particle physics theories, matter at high energies and temperatures 
must be described in terms of helds. Gauge symmetries have proved to be extremely 
useful in describing the standard model of particle physics, according to which at 
high energies the laws of nature are invariant under a nonabelian group G of internal 
symmetry transformations G = SU(3)c x SU(2)l x U{1)y which at a temperature of 
about 200 MeV is spontaneously broken down to G' = SU(3)c x U(l). The subscript 
on the SU(3) subgroup indicates that it is the color symmetry group of the strong 
interactions, SU(2)iX U(l)v is the Glashow-Weinberg-Salam (WS) model of weak 
and electromagnetic interactions, the subscripts L and Y denoting left handedness 
and hypercharge respectively. At low energies, the WS model spontaneously breaks 
to the U(l) subgroup of electromagnetism. 

Spontaneous symmetry breaking is induced by an order parameter (p taking on a 
nontrivial expectation value < p > below a certain temperture T^. In some particle 
physics models, 99 is a fundamental scalar held in a nontrivial representation of the 
gauge group G which is broken. However, p could also be a fermion condensate, as 
in the BGS theory of superconductivity. 
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Earlier we have seen that symmetry breaking phase transitions in gauge held 
theories do not, in general, lead to inhation. In most models, the coupling constants 
which arise in the effective potential for the scalar field (p driving the phase transition 
are too large to generate a period of slow rolling which lasts more than one Hubble 
time Nevertheless, there are interesting remnants of the phase transition: 

topological defects. 

Consider a single component real scalar field with a typical symmetry breaking 
potential 

= (169) 


Unless A -C 1 there will be no inflation. The phase transition will take place on a 
short time scale r < and will lead to correlation regions of radius ^ <t inside of 
which ip is approximately constant, but outside of which p ranges randomly over the 
vacuum manifold Af, the set of values of p which minimizes V{p) - in our example 
p = ± 77 . The correlation regions are separated by domain walls, regions in space 
where p leaves the vacuum manifold M. and where, therefore, potential energy is 
localized. Via the usual gravitational force, this energy density can act as a seed for 
structure. 

Topological defects are familiar from solid state and condensed matter systems. 
Crystal defects, for example, form when water freezes or when a metal crystallizes 0. 
Point defects, linedefects and planar defects are possible. Defects are also common 
in liquid crystals Ea . They arise in a temperature quench from the disordered to the 
ordered phase. Vortices in ^He are analogs of global cosmic strings. Vortices and other 
defects are also produced 0 during a quench below the critical temperature in ^He. 
Finally, vortex lines may play an important role in the theory of superconductivity 0. 

The analogies between defects in particle physics and condensed matter physics 
are quite deep. Defects form for the same reason: the vacuum manifold is topologi¬ 
cally nontrivial. The arguments @ which say that in a theory which admits defects, 
such defects will inevitably form, are applicable both in cosmology and in condensed 
matter physics. Different, however, is the defect dynamics. The motion of defects 
in condensed matter systems is friction-dominated, whereas the defects in cosmology 
obey relativistic equations, second order in time derivatives, since they come from a 
relativistic field theory. 

After these general comments we turn to a classification of topological defects 0 . 
We consider theories with an n-component order parameter p and with a potential 

n 

energy function (free energy density) of the form ( 6 . 1 ) with p'^ = Yl, 


i=l 


There are various types of local and global topological defects (regions of trapped 
energy density) depending on the number n of components of p (see e.g. for a 
comprehensive survey of topological defect models). The more rigorous mathematical 
definition refers to the homotopy of M.. The words “local” and “global” refer to 
whether the symmetry which is broken is a gauge or global symmetry. In the case of 
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local symmetries, the topological defects have a well defined core outside of which (p 
contains no energy density in spite of nonvanishing gradients V(p: the gauge fields 
can absorb the gradient, i.e., D^Lp = 0 when 7 ^ 0 , where the covariant derivative 
is dehned by + ieA^, e being the gauge coupling constant. Global 

topological defects, however, have long range density helds and forces. 

Table 1 contains a list of topological defects with their topological characteristics. 
A “v” marks acceptable theories, a “x” theories which are in conflict with observations 
(for rj ~ 10^® GeV). 


Table 1 

n 

topology 

local defect 

global defect 

domain wall 

1 

no(>M) ^ 1 

X 

X 

cosmic string 

2 

n,(Al) 1 

V 

V 

monopolc 

3 

MM) ^ 1 

X 

V 

texture 

4 

MM) ^ 1 

- 

V 


Theories with domain walls are ruled out since a single domain wall stretching 
across the Universe today would overdose the Universe. Local monopoles are also 
ruled out since they would overdose the Universe. Local textures are ineffective 
at producing structures because there is no traped potential energy. 

From now on we will focus on one type of defects: cosmic strings (see e.g. 
||,[T 0 l,|T 0 l for recent reviews, and |p.04|| for a classic review paper). These arise 
in theories with a complex order parameter (n = 2). In this case the vacuum mani¬ 
fold of the model is 


M = S\ (170) 

which has nonvanishing hrst homotopy group: 

Ui{M) = Z^l. (171) 

A cosmic string is a line of trapped energy density which arises whenever the held (p{x) 
circles Af along a closed path in space {e.g., along a circle). In this case, continuity 
of (p implies that there must be a point with <^9 = 0 on any disk whose boundary is 
the closed path. The points on different sheets connect up to form a line overdensity 
of held energy. 

To construct a held conhguration with a string along the 2 ; axis take 'p{x) to 
cover M. along a circle with radius r about the point {x,y) = ( 0 , 0 ): 
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(^(r, 'd) ~ , r ^ rj 


-1 


(172) 


This configuration has winding number 1, he., it covers Ai exactly once. Maintaining 
cylindrical symmetry, we can extend (|172D to arbitrary r 


Vi.'r, 


(173) 


where /(O) = 0 and /(r) tends to r] for large r. The width w can be found by 
balancing potential and tension energy. The result is 


w ~ 




(174) 


For local cosmic strings, he., strings arising due to the spontaneous breaking of 
a gauge symmetry, the energy density decays exponentially for r S> w. In this case, 
the energy /i per unit length of a string is finite and depends only on the symmetry 
breaking scale rj 


(175) 

(independent of the coupling A). The value of /r is the only free parameter in a cosmic 
string model. 


C. Formation of Defects in Cosmological Phase Transitions 

The symmetry breaking phase transition takes place at T = (called the critical 
temperature). From condensed matter physics it is well known that in many cases 
topological defects form during phase transitions, particularly if the transition rate 
is fast on a scale compared to the system size. When cooling a metal, defects in 
the crystal configuration will be frozen in; during a temperature quench of ^He, thin 
vortex tubes of the normal phase are trapped in the superfiuid; and analogously in a 
temperature quench of a superconductor, flux lines are trapped in a surrounding sea 
of the superconducting Meissner phase. 

In cosmology, the rate at which the phase transition proceeds is given by the 
expansion rate of the Universe. Hence, topological defects will inevitably be produced 
in a cosmological phase transition 0 , provided the underlying particle physics model 
allows such defects. 

The argument which ensures that in theories which admit topological or semitopo- 
logical defects, such defects will be produced during a phase transition in the very 
early Universe is called the Kibble mechanism 0. To illustrate the physics, consider 
a mechanical toy model, first introduced by Mazenko, Unruh and Wald c25l. Take 
(see Fig. 13) a lattice of points on a fiat table. At each point, a pencil is pivoted. 
It is free to rotate and oscillate. The tips of nearest neighbor pencils are connected 
with springs (to mimic the spatial gradient terms in the scalar field Lagrangean). 
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Newtonian gravity creates a potential energy V(ip) for each pencil (cp is the angle 
relative to the vertical direction). V(cp) is minimized for |(p| = rj (in our toy model 
7] = tt/2). Hence, the Lagrangean of this pencil model is analogous to that of a scalar 
held with symmetry breaking potential ( |169|) . 




FIG. 13. The pencil model: the potential energy of a simple pencil has the same form 
as that of scalar fields used for spontaneous symmetry breaking. The springs connecting 
nearest neighbor pencils give rise to contributions to the energy which mimic spatial gradient 
terms in field theory. 


At high temperatures T S> Tc, all pencils undergo large amplitude high frequency 
oscillations. However, by causality, the phases of oscillation of pencils with large 
separation s are uncorrelated. For a system in thermal equilibrium, the length s 
beyond which phases are random is the correlation length ^{t). However, since the 
system is quenched rapidly, there is a causality bound on 

m<t, ( 176 ) 


where t is the causal horizon. 

The critical temperature Tc is the temperature at which the thermal energy is 
equal to the energy a pencil needs to jump from horizontal to vertical position. For 
T < Tc, all pencils want to he flat on the table. However, their orientations are 
random beyond a distance ^{t) determined by equating the free energy gained by 
symmetry breaking (a volume effect) with the gradient energy lost (a surface effect). 
As expected, ^(T) diverges at Tc. Very close to Tc, the thermal energy T is larger 
than the volume energy gain Ecorr in a correlation volume. Hence, these domains 
are unstable to thermal fluctuations. As T decreases, the thermal energy decreases 
more rapidly than Ecorr- Below the “Ginsburg temperature” To, there is insufficient 
thermal energy to excite a correlation volume into the state (p = 0. Domains of size 






(177) 


freeze out olllhl . The boundaries between these domains become topological defects. 
An improved version of this argument has recently been given by Zurek (see also 


Tool ) 
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We conclude that in a theory in which a symmetry breaking phase transitions 
satisfies the topological criteria for the existence of a fixed type of defect, a network of 
such defects will form during the phase transition and will freeze out at the Ginsburg 
temperature. The correlation length is initially given by (|177|) , if the field ip is in 
thermal equilibrium before the transition. Independent of this last assumption, the 
causality bound implies that < ic- 

For times t > to the evolution of the network of defects may be complicated (as 
for cosmic strings) or trivial (as for textures). In any case (see the caveats of ||110|| ), 
the causality bound persists at late times and states that even at late times, the mean 
separation and length scale of defects is bounded by ^{t) < t. 

Applied to cosmic strings, the Kibble mechanism implies that at the time of the 
phase transition, a network of cosmic strings with typical step length ^ (tc) will form. 
According to numerical simulations EdJ , about 80% of the initial energy is in infinite 
strings (strings with curvature radius larger than the Hubble radius) and 20% in 
closed loops. 


D. Evolution of Strings and Scaling 



FIG. 14. Formation of a loop by a self intersection of an infinite string. According to the 
original cosmic string scenario, loops form with a radius R determined by the instantaneous 
coherence length of the infinite string network. 


The evolution of the cosmic string network for t > tc is complicated. The key 
processes are loop production by intersections of infinite strings (see Fig. 14) and 
loop shrinking by gravitational radiation. These two processes combine to create a 
mechanism by which the infinite string network loses energy (and length as measured 
in comoving coordinates). It can be shown (see e.g. |p.04|| ) that, as a consequence, the 
correlation length of the string network is always proportional to its causality limit 


m-t. ( 178 ) 

Hence, the energy density pooit) in long strings is a fixed fraction of the background 
energy density pc(t) 
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p^(t) ~ Pi(t) ^ ~ pt ^ 


( 179 ) 


or 


Poo (t') 
Pc{t) 


~ Gp. 


(180) 


We conclude that the cosmic string network approaches a “scaling solution” in 
which the statistical properties of the network are time independent if all distances 
are scaled to the horizon distance. 


E. Cosmic Strings and Structnre Formation 

The starting point of the structure formation scenario in the cosmic string theory 
is the scaling solution for the cosmic string network, according to which at all times 
t (in particular at teq, the time when perturbations can start to grow) there will be a 
few long strings crossing each Hubble volume, plus a distribution of loops of radius 
R <^t (see Fig. 15). 

The cosmic string model admits three mechanisms for structure formation: loops, 
filaments, and wakes. Cosmic string loops have the same time averaged field as a 
point source with mass M{R) = f3Rp, R being the loop radius and (3 ~ 27r. 
Hence, loops will be seeds for spherical accretion of dust and radiation. 

For loops with R < teq, growth of perturbations in a model dominated by cold 
dark matter starts at teq. Hence, the mass at the present time will be 

M{R, to) = z{teq)(3 Rp . (181) 



FIG. 15. Sketch of the scaling solution for the cosmic string network, 
sponds to one Hubble volume at arbitrary time t. 


The box corre- 
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In the original cosmic string model i’lll it was assnmed that loops dominate over 
wakes. However, according to the newer cosmic string evolntion simnlations IHJ, most 
of the energy in strings is in the long strings, and hence the loop accretion mechanism 
is snbdominant. 

The second mechanism involves long strings moving with relativistic sp eed in 
their normal plane which give rise to velocity pertnrbations in their wake The 
mechanism is illnstrated in Fig. 16: space normal to the string is a cone with dehcit 
angle 0 


a = SirGjj,. 


( 182 ) 


If the string is moving with normal velocity v throngh a bath of dark matter, a 
velocity pertnrbation 


Sv = A'nGfxv'y 


( 183 ) 


[with 7 = (1 — towards the plane behind the string resnlts. At times af¬ 

ter teq, this indnces planar overdensities, the most prominent (he., thickest at the 
present time) and nnmerons of which were created at t^q, the time of eqnal mat¬ 
ter and radiation EiZrtiJ. The corresponding planar dimensions are (in comoving 
coordinates) 


teqZ{teq) X teqZ(teq)V ~ (40 X 40t>) MpC^ . 


(184) 



I 

FIG. 16. Sketch of the mechanism by which a long straight cosmic string S moving 
with velocity v in transverse direction through a plasma induces a velocity perturbation At> 
towards the wake. Shown on the left is the deficit angle, in the center is a sketch of the 
string moving in the plasma, and on the right is the sketch of how the plasma moves in the 
frame in which the string is at rest. 


The thickness d of these wakes can be calcnlated using the Zel’dovich 
approximation The result is (for G/i = 10“®) 

d ~ Gjjiv^{y)z{teqY teq — 4"^ Mpc . (185) 

Wakes arise if there is little small scale structure on the string. In this case, the 
string tension equals the mass density, the string moves at relativistic speeds, and 
there is no local gravitational attraction towards the string. 
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In contrast, if there is small scale structure on strings, then the string tension T is 
smaller 13 than the mass per unit length /r , and thus there is a gravitational force 
towards the string which gives rise to cylindrical accretion, producing hlaments 
Which of the mechanisms - hlaments or wakes - dominates is determined by the 
competition between the velocity induced by the Newtonian gravitational potential 
of the strings and the velocity perturbation of the wake. 

The cosmic string model predicts a scale-invariant spectrum of density pertur¬ 
bations, exactly like inhationary Universe models but for a rather diherent reason. 
Consider the r.m.s. mass huctuations on a length scale 2'Kk~^ at the time tuik) when 
this scale enters the Hubble radius. From the cosmic string scaling solution it follows 
that a hxed (he., tuik) independent) number v of strings of length of the order tnik) 
contribute to the mass excess SM{k, tnik)). Thus 


6M 


{k, tnik)) 


vntnik) 


G-Hj,^{k)tUk) 


V Gfx. 


(186) 


Note that the above argument predicting a scale invariant spectrum will hold for all 
topological defect models which have a scaling solution, in particular also for global 
monopoles and textures. 

The amplitude of the r.m.s. mass huctuations (equivalently: of the power spec¬ 
trum) can be used to normalize Gfi. Since today on galaxy cluster scales 


6M 

W 


{k, to) 


rsj 


1 , 


the growth rate of huctuations linear in a{t) yields 


6M 


{k, t, 


eqj 


10 


-4 


and therefore, using h ~ 10, 




(187) 


(188) 


(189) 


A similar value is obtained by normalizing the model to the COBE amplitude of CMB 
anisotropies on large angular scales (the normalizations from COBE and from 

the power spectrum of density perturbations on large scales agree to within a factor 
of 2). Thus, if cosmic strings are to be relevant for structure formation, they must 
arise due to a symmetry breaking at energy scale r] ~ lO^^GeV. This scale happens to 
be the scale of unihcation (GUT) of weak, strong and electromagnetic interactions. 
It is tantalizing to speculate that cosmology is telling us that there indeed was new 
physics at the GUT scale. 

A big advantage of the cosmic string model over inhationary Universe models 
is that HDM is a viable dark matter candidate. Gosmic string loops survive free 
streaming, as discussed in Section 3.B, and can generate nonlinear structures on 
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galactic scales, as discussed in detail in ||125|,|126|| . Accretion of hot dark matter by a 
string wake was studied in [Id^]. In this case, nonlinear perturbations develop only 
late. At some time tnu all scales up to a distance g^ax from the wake center go 
nonlinear. Here 


gmax ~ G^V-f{v)z{teg)Heq ~ 4n MpC , (190) 

and it is the comoving thickness of the wake at tni. Demanding that tni corresponds 
to a redshift greater than 1 leads to the constraint 

G'/r>5■10■^ (191) 


Note that in a cosmic string and hot dark matter model, wakes form nonlinear struc¬ 
tures only very recently. Accretion onto loops and small scale structure on the long 
strings provide two mechanisms which may lead to high redshift objects such as 
quasars and high redshift galaxies. The hrst mechanism has recently been studied in 
1^1, the second in ||128| , |129|| . 


The power spectrum of density fluctuations in a cosmic string model with HDM 
has recently been studied numerically by Mahonen based on previous work of 
El (see also ||132|| for an earlier semi-analytical study). The spectral shape agrees 
quite well with observations, and a bias factor of less than 2 is required to give the 
best-£t amplitude for a CODE normalized model. Note, however, that the results 
depend quite sensitively on the details of the string scaling solution which are at 
present not well understood. 

Due to lack of space, I will not discuss the global monopole [^231 and global 
texture models of structure formation. There has been a lot of work on the 
texture model, and the reader is referred to ||136|,|135|| for recent review articles. 


F. Signatures for Strings 


The cosmic string theory of structure formation makes several distinctive predic¬ 
tions, both in terms of the galaxy distribution and in terms of CMB anisotropies. On 
large scales (corresponding to the comoving Hubble radius at teq and lar ger , struc¬ 


ture is predicted to be dominated either by planar lliZrEdi or filamentary galaxy 
concentrations. For models in which the strings have no local gravity, the resulting 
nonlinear structures will look very different from the nonlinear structures in models 
in which local gravity is the dominant force. As discovered and discussed recently in 
a baryon number excess is predicted in the nonlinear wakes. This may explain 
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the “cluster baryon crisis” 123, the fact that the ratio of baryons to dark matter in 
rich clusters is larger than what is compatible with the nucleosynthesis constraints in 
a spatially flat Universe. 
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FIG. 17. Sketch of the Kaiser-Stebbins effect by which cosmic strings produce linear 
discontinuities in the CMB. Photons 7 passing on different sides of a moving string S 
(velocity v) towards the observer O receive a relative Doppler shift due to the conical 
nature of space perpendicular to the string (deficit angle a). 


As described in the previous subsection, space perpendicular to a long straight 
cosmic string is conical with deficit angle given by (p.82|) . Consider now CMB radiation 
approaching an observer in a direction normal to the plane spanned by the string and 
its velocity vector (see Fig. 17). Photons arriving at the observer having passed on 
different sides of the string will obtain a relative Doppler shift which translates into 
a temperature discontinuity of amplitude 

6T 

— = AnGuv-fiv ), (192) 


where v is the velocity of the string. Thus, the distinctive signature for cosmic strings 
in the microwave sky are line discontinuities in T of the above magnitude. 

Given ideal maps of the CMB sky it would be easy to detect strings. However, 
real experiments have hnite beam width. Taking into account averaging over a scale 
corresponding to the beam width will smear out the discontinuities, and it turns out 
to be surprisingly hard to distinguish the predictions of the cosmic string model from 
that of inflation-based theories using quantitative statistics which are easy to evaluate 
analytically, such as the kurtosis of the spatial gradient map of the CMB . There 
may be ways to distinguish between string and inflationary models by looking at the 
angular power s p ectr um of CMB anisotropies. Work on this subject, however, is still 
controversial fiHiH. 

Global textures also produce distinctive non-Gaussian signatures Edil in CMB 
maps. In fact, these signatures are more pronounced and on larger scales than the 
signatures in the cosmic string model. 
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G. Principles of Baryogenesis 


Baryogenesis is another area where particle physics and cosmology connect in a 
very deep way. The goal is to explain the observed asymmetry between matter and 
antimatter in the Universe. In particnlar, the objective is to be able to explain the 
observed valne of the net baryon to entropy ratio at the present time 

^(to) ~ 10-^° (193) 

s 

starting from initial conditions in the very early Universe when this ratio vanishes. 
Here, Ahb is the net baryon nnmbfir density and s the entropy density. 

As pointed ont by Sakharov , three basic criteria mnst be satished in order to 
have a chance at explaining the data: 


1. The theory describing the microphysics mnst contain baryon nnmber violating 
processes. 

2 . These processes mnst be C and CP violating. 


3. The baryon nnmber violating processes mnst occnr ont of thermal eqnilibrinm. 


As was discovered in the 1970’s Eda, all three criteria can be satished in GUT 
theories. In these models, baryon nnmber violating processes are mediated by sn- 
perheavy Higgs and gange particles. The baryon nnmber violation is visible in the 
Lagrangian, and occnrs in pertnrbation theory (and is therefore in principle easy to 
calcnlate). In addition to standard model CP violation, there are typically many 
new sonrces of CP violation in the GUT sector. The third Sakharov condition can 
also be realized: After the GUT symmetry-breaking phase transition, the snperheavy 
particles may fall ont of thermal eqnilibrinm. The ont-of-eqnilibrinm decay of these 
particles can thns generate a nonvanishing baryon to entropy ratio. 

The magnitnde of the predicted ub/s depends on the asymmetry e per decay, on 
the conpling constant A of the ub violating processes, and on the ratio nx/s of the 
nnmber density nx of snperheavy Higgs and gange particles to the nnmber density 
of photons, evalnated at the time td when the baryon nnmber violating processes 
fall ont of thermal eqnilibrinm, and assnming that this time occnrs after the phase 
transition. The qnantity e: is proportional to the CP-violation parameter in the model. 
In a GUT theory, this GP violation parameter can be large (order 1), whereas in the 
standard electroweak theory it is given by the GP violating phases in the GKM mass 
matrix and is very small. As shown in ||146|| it is easily possible to constrnct models 
which give the right ub/ s ratio after the GUT phase transition (for recent reviews of 
baryogenesis see ||147|| and |p.48|| ). 
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H. GUT Baryogenesis and Topological Defects 

The ratio ub/s, however, does not only depend on e, but also on nx/s(td). If the 
temperature T^, at the time td is greater than the mass nix of the superheavy particles, 
then it follows from the thermal history in standard cosmology that nx ~ s. However, 
if Td < mx, then the number density of X particles is diluted exponentially in the 
time interval between when T = mx and when T = T^. Thus, the predicted baryon 
to entropy ratio is exponentially suppressed: 

— ~ (194) 

s g* 

where g* is the number of spin degrees of freedom in thermal equilibrium at the time 
of the phase transition. In this case, the standard GUT baryogenesis mechanism is 
ineffective. 

However, topological defects may come to the rescue ^0. As we have seen ear¬ 
lier in this section, topological defects will inevitably be produced in the symmetry 
breaking GUT transition provided they are topologically allowed in that symmetry 
breaking scheme. The topological defects provide an alternative mechanism of GUT 
baryogenesis. 

Inside of topological defects, the GUT symmetry is restored. In fact, the defects 
can be viewed as solitonic conhgurations of X particles. The continuous decay of 
defects at times after td provides an alternative way to generate a nonvanishing baryon 
to entropy ratio. The defects constitute out of equilibrium conhgurations, and hence 
their decay can produce a nonvanishing ub/s in the same way as the decay of free X 
quanta. 

The way to compute the resulting tib/s ratio is as follows: The defect scaling 
solution gives the energy density in defects at all times. Taking the time derivative 
of this density, and taking into account the expansion of the Universe, we obtain the 
loss of energy attributed to defect decay. By energetics, we can estimate the number 
of decays of individual quanta which the defect decay corresponds to. We can then 
use the usual perturbative results to compute the resulting net baryon number. 

Provided that mx < Td, then at the time when the baryon number violating 
processes fall out of equilibrium (when we start generating a nonvanishing hb) the 
energy density in free X quanta is much larger than the defect density, and hence 
the defect-driven baryogenesis mechanism is subdominant. However, if mx > Td, 
then as indicated in (|194| ), the energy density in free quanta decays exponentially. In 
constrast, the density in defects only decreases as a power of time, and hence soon 
dominates baryogenesis. 

One of the most important ingredients in the calculation is the time dependence 
of ^{t), the separation between defects. Immediately after the phase transition at the 
time f/ of the formation of the defect network, the separation is ^(t/) ~ In 

the time period immediately following, the time period of relevance for baryogenesis, 
^{t) approaches the Hubble radius according to the equation 
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( 195 ) 


e(() c. e((,)(f )=-'*. 

Using this result to calculate the defect density, we obtain after some algebra 


ns, ,2 I 

defect 0 ? 

S T] S 


(196) 


where ub/sIo is the unsuppressed value of ub/s which can be obtained using the 
standard GUT baryogenesis mechanism. We see from (1^) that even for low values 
of Td, the magnitude of ub/s which is obtained via the defect mechanism is only 
suppressed by a power of T^,. However, the maximum strength of the defect channel 
is smaller than the maximum strength of the usual mechanism by a geometrical sup¬ 
pression factor which expresses the fact that even at the time of defect formation, 
the defect network only occupies a small volume. 


I. Electroweak Baryogenesis and Topological Defects 


It has been known for some time that there are baryon number violating processes 
even in the standard electroweak theory. These processes are, however, nonpertur- 
bative. They are connected with the t’Hooft anomaly which in turn is due to 
the fact that the gauge theory vacuum is degenerate, and that the different degen¬ 
erate vacuum states have different quantum numbers (Chern-Simons numbers). In 
theories with fermions, this implies different baryon number. Conhgurations such as 
sphalerons I® which interpolate between two such vacuum states thus correspond to 
baryon number violating processes. 

As pointed out in | 152 |, the anomalous baryon number violating processes are in 
thermal equilibrium above the electroweak symmetry breaking scale. Therefore, any 
net baryon to entropy ratio generated at a higher scale will be erased, unless this 
ratio is protected by an additional quantum number such as a nonvanishing B — L 
which is conserved by electroweak processes. 

However, as hrst suggested in ||153|| and discussed in detail in many recent papers 


(see |p.54|| for reviews of the literature), it is possible to regenerate a nonvanishing 
hb/s below the electroweak symmetry breaking scale. Since there are hb violating 
processes and both C and CP violation in the standard model, Sakharov’s conditions 
are satished provided that one can realize an out-of-equilibrium state after the phase 
transition. Standard model CP violation is extremely weak. Thus, it appears neces¬ 
sary to add some sector with extra CP violation to the standard model in order to 
obtain an appreciable ub/s ratio. A simple possibility which has been invoked often 
is to add a second Higgs doublet to the theory, with CP violating relative phases. 

The standard way to obtain out-of-equilibrium baryon number violating processes 
immediately after the electroweak phase transition is to assume that the transi¬ 
tion is strongly hrst order and proceeds by the nucleation of bubbles (note that these 
are two assumptions, the second being stronger than the hrst!). 
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Bubbles are out-of-equilibrium configurations. Outside of the bubble (in the false 
vacuum), the baryon number violating processes are unsuppressed, inside they are 
exponentially suppressed. In the bubble wall, the Higgs helds have a nontrivial pro- 
hle, and hence (in models with additional CP violation in the Higgs sector) there is 
enhanced CP violation in the bubble wall. In order to obtain net baryon production, 
one may either use fermion scattering off bubble walls EH (because of the CP viola¬ 
tion in the scattering, this generates a lepton asymmetry outside the bubble which 
converts via sphalerons to a baryon asymmetry) or sphaleron processes in the bubble 


wall itself . It has been shown that, using optimistic parameters (in particular 


a large CP violating phase A9cp in the Higgs sector) it is possible to generate the 
observed np/s ratio. The resulting baryon to entropy ratio is of the order 


Hb 


~ a 


w 


(9^ 




Aft 


CP 
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where aw refers to the electroweak interaction strength, g* is the number of spin 
degrees of freedom in thermal quilibrium at the time of the phase transition, and 
is the top quark mass. The dependence on the top quark mass enters because net 
baryogenesis only appears at the one-loop level. 

However, analytical and numerical studies show that, for the large Higgs masses 
which are indicated by the current experimental bounds, the electroweak phase tran¬ 
sition will unlikely be sufficiently strongly hrst order to proceed by bubble nucleation. 
In addition, there are some concerns as to whether it will proceed by bubble nucleation 
at all (see e.g. ||15^| ). 

Once again, topological defects come to the rescue. In models which admit de¬ 
fects, such defects will inevitably be produced in a phase transition independent of 
its order. Moving topological defects can play the same role in baryogenesis as nucle¬ 
ating bubbles. In the defect core, the electroweak symmetry is unbroken and hence 
sphaleron processes are unsuppressed In the defect walls there is enhanced CP 
violation for the same reason as in bubble walls. Hence, at a fixed point in space, a 
nonvanishing baryon number will be produced when a topological defect passes by. 
Defect-mediated electroweak baryogenesis has been worked out in detail in ||160| 


(see for previous work) in the case of cosmic strings. The scenario is as follows: 
at a particular point x in space, antibaryons are produced when the front side of the 
defect passes by. While x is in the defect core, partial equilibration of takes place 
via sphaleron processes. As the back side of the defect passes by, the same number of 
baryons are produced as the number of antibaryons when the front side of the defect 
passes by. Thus, at the end a positive number of baryons are left behind. 

As in the case of defect-mediated GUT baryogenesis, the strength of defect- 
mediated electroweak baryogenesis is suppressed by the ratio SF of the volume which 
is passed by defects divided by the total volume, i.e. 


Ub 

s 


SF^Io 

s 


(198) 
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where {ub/s)\q is the result of (|197|) obtained in the bubble nucleation mechanism. 

A big caveat for defect-mediated electroweak baryogenesis is that the standard 
electroweak theory does not admit topological defects. However, in a theory with 
additional physics just above the electroweak scale it is possible to obtain defects 
(see e.g. ||162|| for some specihc models). The closer the scale r] of the new physics is 
to the electroweak scale rjEw, the larger the volume in defects and the more efficient 
defect-mediated electroweak baryogenesis. Using the result of ( |195D for the separation 
of defects, we obtain 


SF ~ (199) 

V 

Obviously, the advantage of the defect-mediated baryongenesis scenario is that it 
does not depend on the order and on the detailed dynamics of the electroweak phase 
transition. 


J. Summary 


As we have seen, topological defects may play an important role in cosmology. 
Defects are inevitably produced during symmetry breaking phase transitions in the 
early Universe in all theories in which defects are topologically stable. Theories giving 
rise to domain walls or local monopoles are ruled out by cosmological constraints. 
Those producing cosmic strings, global monopoles and textures are quite attractive. 

If the scale of symmetry breaking at which the defects are produced is about 
10^® GeV, then defects can act as the seeds for galaxy formation. Defect theories of 
structure formation predict a roughly scale-invariant spectrum of density perturba¬ 
tions, similar to inflation-based models. However, the phases in the density held are 
distributed in a non-Gaussian manner. Thus, the predictions of defect models can 
be distinguished from those of inhationary models. In addition, the predictions of 
diherent defect models can be distinguished from eachother. 

Focusing on the cosmic string theory of structure formation, we have seen that 
the model gives rise to several distinctive signatures. The large-scale structure in the 
Universe is predicted to be dominated by either planar or hlamentary structures (de¬ 
pending on whether there is small-scale structure on the strings), with a distinctive 
scale given by the comoving Hubble radius at teq- There are many ways topolog¬ 
ical measures with which one can quantify this prediction (see e.g. ||163|| ). In the 
GMB, the distinctive signature for strings is predicted to be line discontinuities in the 
temperature maps. 

Topological defects may also play a crucial role in baryogenesis. This applies both 
to GUT and electroweak baryogenesis. The crucial point is that defects constitute out- 
of-equilibrium conhgurations, and may therefore be the sites of net baryon production. 
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